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ABSTRACT 

Using the H 13 CN and HN 13 C J=l-0 line observations, the abundance ratio 
of HCN/HNC has been estimated for different evolutionary stages of massive 
star formation: Infrared dark clouds (IRDCs), High-mass protostellar object 
(HMPOs), and Ultra-compact HII regions (UCHIIs). IRDCs were divided into 
‘quiescent IRDC cores’ and ‘active IRDC cores’, depending on star formation 
activity. The HCN/HNC ratio is known to be higher at active and high temper¬ 
ature regions related to ongoing star formation, compared to cold and quiescent 
regions. Our observations toward 8 quiescent IRDC cores, 16 active IRDC cores, 
23 HMPOs, and 31 UCHIIs show consistent results; the ratio is 0.97 (± 0.10), 
2.65 (± 0.88), 4.17 (± 1.03) and 8.96 (± 3.32) in these respective evolution¬ 
ary stages, increasing from quiescent IRDC cores to UCHIIs. The change of 
the HCN/HNC abundance ratio, therefore, seems directly associated with the 
evolutionary stages of star formation, which have different temperatures. One 
suggested explanation for this trend is the conversion of HNC to HCN, which 
occurs effectively at higher temperatures. To test the explanation, we performed 
a simple chemical model calculation. In order to fit the observed results, the 
energy barrier of the conversion must be much lower than the value provided by 
theoretical calculations. 


1. Introduction 

Massive stars are believed to have a decisive effect on the physical and chemical 
evolution of galaxies through their energetic feedbacks, including ultraviolet (UV) 
radiation, stellar winds, bipolar outflows, HII regions, and eventually supernova explosions. 
Nevertheless, our understanding about how massive stars form is not yet well established (see 
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Zinnc cker fe York 2007 for a review). It is hard to capture the early evolutionary phases 
of massive stars due to their fast evolution. In addition, observational information can be 
easily contaminated because the objects are distant and typically form in clusters. For 
these reasons, observational and theoretical works regarding massive star formation are 
still relatively limited. Therefore, characterizing the various early phases of massive star 
formation in a coherent way is important to better understand massive star formation. 

HCN and its geometrical isomer HNC are two fundamental and ubiquitous molecules 


in the dense interstellar medium. They appear to be a tracer of high-c 


HCN has b een used as a tracer of infall motion in star forming regions flWu & Ev ans 


Sqhnetah 


ensity gas, and 


2003 


2007). Above all, since t 


depend on the kinetic temperature 


he HCN/HNC at 

mndance ratio is known 

to strongly 

), it can 

Goldsmith et al. 

1986; 

Schilke et al. 

1992 


be a useful tool to explore the physical and chemical conditions of star-forming regions. 
Although both HCN and HN C are k nown to form in equal measure through the following 


dissociative recombination (Mendes et al 


2012 ), 


HCNH+ + e 


HCN + H 


( 1 ) 


HNC + H 

the HCN/HNC abundance ratio has been shown to vary by regions. For example, 
observational studies of the Orion Molecular Cloud (OMC-1) sh owed that the abundance 


ratio increases from 5 at the cold edge to 80 in the warm core flSchilke et al. 


19921 ). To 


explain such a variation in the HCN/HNC abundance ratio, many mechanisms have been 
suggested that the neutral-neutral reaction below, which selectively consumes HNC, is 
likely the main cause of the temperature dependence of the HCN/HNC abundance ratio. 

HNC + H ->• HCN + H. (2) 


Schilke et al. 


(119921 1 hypothesized that the above reaction is only effective at high 


temperature, and empirically set the rate coefficient to reproduce the variance of the 
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HCN/HNC abundance ratio in OMC-1. In addition, iTalbi et al. f 19961 ) theoretically 
analyzed the neutral-neutral reaction using quantum-chemical calculation, and the 


resulting rate coefficient for the neutral-n eutral reaction has 


UMIST ((McElroy ct ah 


2013 ) and KID A (Harada et al. 


2010 , 


networks with the associated conventional chemical models. 


leen reflected in the 


2012) chemical reaction 


Talbi et al. 


01996! ) suggested 


that an activation energy barrier for the reaction is as high as 2000 K so that the reaction 
is not effective at a temperature lower than 100 K. However, this explanation does not fully 


account for t 


lmhuctuation of theJdCN/HNC abundance ratio in observations. For example, 


Schilke et al 


dataset of 


(1992) and 


Schilke et al. 


Hirota et al. 


( 1998 ) fitted the HCN/HNC abundance ratio for the 


(119921 ) by modifying the rate coefficient in a semi-quantitative way. 


Based on the fitting, they found that the temperature dependence of the abundance ratio 
still holds at T > 24K, and the rate coefficient, for which the activation energy is set to 
be ~ 200 K, could well explain the observational results. Such a discre pancy between th e 


obser vational works and the theoretical works still remains unresolved (jGraninger et al. 


2014 ). 


There have been several attem pts to measure t 


massive star forming regions (IVasvunina et al. 


2011 


tie HCN/H 


Hoq et al 


C abundance ratio toward 


2013 : 


Gerner et al. 


20141). 


However, many of those studies dealt with limited phases of massive star formation or 
suffered from the optical depth effect in the calculation of abundances. In this study, we 
determine the HCN/HNC abundance ratio based on optically thin lines toward various 
evolutionary stages related to massive star formation; infrared dark clouds (IRDCs), 
high-mass protostellar objects (HMPOs) and ultracompact HII regions (UCHIIs). 

IRDCs were identified by the MSX and Spitzer s urveys as dark extinction features 


against the bright Ga 


Peretto fe Fuller 


actic mid-infrared background (lEgan et al. 


1998 : 


Simon et al. 


2006a ; 


20091 ). Since these objects contain cold (<25 K) and dense (> 10 5 cm 


- 3 \ 
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cores with strong (sub)millimeter emissions (e.g., Rathborne et al. 2006), they are thought 
to be the ideal birthplaces of massive stars. HMPOs are luminous infrared (Lhoi > 10 3 


point-li ke sources wit ho ut associatec 


Sridharan et al. 


2002 


Beuther et ah 


radio continuum emission flMolinari et al. 


1996 


2000 


2002). They have highe r gas temperatures than IRDCs 


and mostly show outflow activities (e.g., 


(n P > 10_cm 3 ), 


Zhang et al. 


( 20051 )). Thus they are believed to be 


1989; 


Kurtz et al. 


1994 


Garcia-Segura et al 


1996 


ionized regions ( 

Wood & Churchwell 

Kim & Koo 

2001 

). The objects are 


considered to represent the childhood of HII regions. They are the most evolved among the 
three phases discussed here but are still relatively early phases of massive star formation. 


This paper is organized as follows. The details about source selection, observations, 
and used archive data are presented in § 2. The observed spectra and the derivation of 
column densities of HCN and HNC and the abundance ratio between them are provided 
in § 3. The analysis for the derived HCN/HNC abundance ratio and associated chemical 
modeling are discussed in § 4. The main results are summarized in § 5. 


Observations and archival data 


2.1. Target selection 


After 38 IRDCs were mapped in 1 . 2 mm continuum emiss i on an d identified into 


190 compact cores by 


Rathborne et al. 


(120061) . 


Chambers et al 


(120091 ) classified them as 


‘quiescent’ prestcllar cores and ‘active’ protostellar cores by examining the star formation 
activity in the cores. Active IRDC cores (hereafter aIRDCc) contain both 4.5 and 24 /an 
infrared emissions, known as the signature of star forming activities, while quiescent IRDC 
cores (hereafter qIRDCc) show neither emission. Hereafter, those IRDC cores are treated 
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in analysis as quiescent or active because each sub-sample shows these separate properties 
in our observati onal results. We adopted 19 qIRDCc and 35 aIRDCc from the catalog of 


Chambers et al. 


020091 ) as our IRDC targets with following conditions: no maser sources, 


no nearby (32") infrared sources, and well isolated from other adjacent cores. However, 
because of a low detection rate for the qIRDCc, it was necessary to increase the number of 
quiescent objects, and 13 qIRDCc were added. For the additional targets, we preferentially 
selected massive cores (> 70M 0 ) because qIRDCc detected both in the H 13 CN and HN 13 C 
lines had their mass greater than ~ 1OOM 0 . 


We selected 69 HMPOs from the catalogs of 


( 1996 ) and 54 UCHIIs from the catalogs of Wood fe C hurchwelll ( 19891) and 


Sridharan et al. 


( 2002 ) and 


Molinari et al. 


Kurtz et ah 


(119941) . These selected targets have their (sub)millimeter continuum mapping data, which 


can be used for the H 2 column density calculation. For 


MAMBO 1.2 mm continuum maps from 


Beuther et al. 


HMPOs, we were able to use the 


(120021 ). Consequently, our sample 


consisted of 67 IRDCs (32 qIRDCc and 35 aIRDCc), 69 HMPOs, and 54 UCHIIs. 


2.2. Observations 


The J=l—0 transitions of HCN, H 13 CN, HNC, and HN 13 C were observed in 2012 
May, 2012 November, 2013 February, and 2 014 Jan uar y and Mar ch using th e KVN Yonsei 


2011 


Lee et al. 


20111) . The lines 


21 m and the KVN Ulsan 21 m telescope (Kim et al. 
were observed with the KVN 86 GHz SIS mixer receiver, which covers a frequency range of 
(85—95) GHz. The rest frequencies, dipole moments, and intrinsic line strengths of observed 
lines are summarized in Table CD The velocity resolution of all lines is 0.43 km s _1 . All 
observations were carried out with the position switching mode. The on-source integration 
time of the HCN and HNC lines for all sources was 10 minutes, resulting in a rms noise 
tempertature of ~50 mK. For the isotopic lines, all but 13 additional qIRDCc had on-source 
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integration time of 30 minutes. For the 13 qIRDCc, we gave an on-source integration time of 
~1 hour because these cores appear faint in H 13 CN compared to other sub-samples. These 
sources were only observed in H 13 CN and HN 13 C J = 1-0. The main-beam efficiencies are 
0.43 and 0.37 for the KVN Yonsei and Ulsan telescopes, respectively, and the half-power 
beam width (HPBW) of both telescopes is 32". The system temperature ranges from 170 K 
to 280 K, and the pointing and focus were adjusted by observing strong SiO maser sources 
every one to two hours. The line intensities were calibrated by the standard chopper wheel 
method. Except for the 13 qIRDCc observation, we first observed all the sources in HCN 
and HNC J=1-0, then observed H 13 CN and HN 13 C J=l-0 lines only for the sources with 
-Z~a*, P eak > ~0.5 K in HCN and HNC J=l-0. All spectra were reduced using CLASS in the 
GILDAS software package. 


2.3. Archival data 


In order to determine the fractional abundance of HCN(or HNC) relative to H 2 , the 


H 2 column density must 


map data flBeuther et al. 


je ob tained. For this, we used the MAMBO 1.2 mm continuum 


SCUBA legacy catalogue, 


2002) and the SCUBA 85 0 /im continuum map data (the JCMT 


Di Francesco et al. 


(120081 ) j. 


The MAMBO (37-element array) 1.2 mm data were obtained with the dua . 


mapping mode and the maps typically have a resolution of 11" flBeuther et al 


-beam OTF 


2002) . The 


SCUBA Legacy Catalogues consist of data which exists in the JCMT archive. Because a 
large number of IRDCs, HMPOs, and UCHIIs in our sample have corresponding SCUBA 
objects counterparts, we can use these homogenous datasets throughout all three types of 
objects. The HPBW of JCMT is 13.5" at 850 pm. However, the SCUBA Legacy Catalogue 
data have a larger effective HPBW of 22.9"; this is not only because of the effective 
smoothing but also because of the additional smoothing applied to reduce pixel-to-pixel 





















noise |Di Francesco et ah 


20081 ). The flux uncertainties of both map data sets are ~20 %. 


3. Results and analysis 
3.1. Observed spectra 


We detected both H 13 CN and HN 13 C lines toward 8 qIRDCc, 16 aIRDCc, 23 HMPOs 
and 31 UCHIIs with 3-cr criterion. The information for these sources are listed in Tables [2H5J 
For the H 13 CN line, the sources with the F = 2-1 hyperfine component intensity greater 
than the signal-to-noise ratio (SNR) of 3 are regarded to be detected. However, because of 
the simple line feature of HN 13 C, the detection of the HN 13 C line in three sources could be 
confirmed by eye even though their SNR is smaller than 3. Tables [6H9] present the observed 
line parameters for qIRDCc, aIRDCc, HMPOs, and UCHIIs, respectively. All the line 
parameters of H 13 CN are derived from the multiple Gaussian fitting while those of HN 13 C 
are determined with a single Gaussian fit. Figure |T] exhibits the representative spectra of 
each sub-sample. The HCN spectra show three hyperfine components originating from the 
nuclear spin interaction. In optically thin and LTE conditions, the relative intensities of 
F = 2-1, F = 1-1, F = 1-0 hyperfine transitions of HCN should be 5 : 3 : 1, reflecting 
their statistical weights. However, our HCN and HNC lines seem to be optically thick and 
highly blended. Also, the HCN and H 13 CN lines show an ano m alous hyperfi ne s tructure 


from what would be expected under the LTE condition (ILoughnane et al. 


20121). The 


HNC and HN 1 3 C li nes also h ave a hyp erfine structure, but their splittings are too small 


^0.7 km s 


-l. 


van der Tak et al. 


( 20091 )) to be resolved in our sources. 


Many studies have reported such HCN anomalies and 

lave struggled to examine 

the underlying mechanism for those HCN anomalies fe.g.. 

Cernicha 

iro et al. 

1984 


Cernicharo & Guelin 

1987; 

Louffhnane et al. 

2012 

)• 

Louehnane et al. 

(2012 

) concluded that 
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HCN hyperfine anomalies are common in both low- and high-mass star-forming regions 
and inferred line overlap effects as the main origins of the anomalies. They have pointed out 
that the anomalies possibly lead to significant errors in estimating opacities with hyperfine 
structure fitting. As a result, the H 13 CN and HN 13 C lines were only chosen for quantitative 
measure of column density under the assumption that both lines are optically thin. 

Figure [2] shows the averaged spectra of H 13 CN and HN 13 C for each sub-sample. The 
peak intensity of the averaged H 13 CN spectrum increases from qIRDCc to UCHIIs, while 
the HN 13 C intensity remains approximately constant, indicative of the increase of the 
HCN/HNC abundance ratio with the evolutionary stage. The averaged spectra of aIRDCc 
and HMPOs appear to have similar peak intensities and line widths. The detailed analysis 
of the HCN/HNC abundance ratio will be discussed in §4.1. 


3.2. Column densities of HCN, HNC, and H 2 


To determine the HCN/HNC abundance ratio, we assume that the observed H 13 CN J 
= 1—0 and HN 13 C J = 1—0 lines arise from the same region. In ord er to test the reliability 


of the assumption, we compare the widths of two lines (ISakai et al. 


20101 ): the HN 13 C J = 


1-0 line is fitted with a single Gaussian profile while the H 13 CN J = 1-0 line is fitted with 
its hyperfine structure to calculate the line width. In some cases, the hyperfine structure 
fit to the H 13 CN line was not reliable because the transitions are highly blended, and such 
sources were excluded by eye from the line width comparison. 


Figure [3] shows the derived line width of HN 13 C J = 1—0 against that of H 13 CN 
J = 1—0. All evolutionary phases show a correlation between two line widths in the 
confidence level above 99 %. These correlations indicate that both lines arise from the same 
region within a clump. In this figure, we plot qIRDCc and aIRDCc together because the 
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number of qIRDCc sources that can be used for this correlation is too small to provide a 
reasonable test. Above all, according to various mapping surveys of massive star forming 


regions including qIRDCc, the HCN and h 


targets flBergin et al. 


19971 : 


Hoq et al. 


2013 


MC lines show similar distributions in many 


Miettinen et ah 


2014 : 


Rathborne et al. 


2014!). 


This consistent spatial distribution of two line emission supports that our assumption is 
reasonable and thus, the HCN/HNC abundance ratio can be used to study the physical 
properties of each source. 


We also assume that the observed H 13 CN J = 1—0 and HN 13 C J = 1—0 lines are 
optically thin. In order to test the reliability of the assumption, the optical depth of 
these lines should be obtained. However, as mentioned in §3.1, opacities derived from the 
hyperhne structure fitting are not reliable due to the anomalous H 13 CN hyperhne structure. 
In addition, the hyperhne structure of HN 13 C J=l-0 is not resolved in our sources. As 
a result, we derived the optical depth and the excitation temperature by applying the 
standard LTE approximation for simplicity. 


If the optically thick main isotopic line is observed, the optical depth of the H 13 CN line 
can be derived from 


Ahin = - In 


T* 


thin 


( 3 ) 


[T ex — J u (T hg )\ 

where T t hm is the beam corrected brightness temperature, T ex is the excitation temperature 
derived from the optically thick main isotopic line, and Tb g is the background radiation 
temperature. Since the HCN J = 1-0 emission is optically thick, we can calculate an 
excitation temperature according to 


T = 

fiX 


hv u 


k 


In 1 + 


( hvjk ) 


Ti 


thick 


Jv(T' bg) 


-1 


( 4 ) 


where Tthick is the beam corrected brightness temperature of the HCN line fjPurcell et al. 


20061 ). The same calculation has been applied to the HNC and HN 13 C lines. The derived 
values for each source are listed in Tables [6] [9] Also, the mean values and standard 
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deviations of r t hi n and T ex derived from individual evolutionary stages are summarized in 
Table [TUI According to the results, the r t hi n are as low as ~ 0.1, and thus the optically thin 
assumption for both isotopologues seems to be reasonable. 

The derived T ex is considerably lower than the temperature quoted in previous 
work. Such low T ex may suggest that the rotational levels of these molecules are excited 
sub-thermally. For the total column density of H 13 CN or HN 13 C, we adopt these excitation 
temperatures derived from the optically thick HCN or HNC lines as our standard; all level 
populations of a given molecule are determined by one excitation temperature although 
these molecules are sub-thermally excited. The detailed non-LTE analysis is beyond 
the scope of this work. However, we note that our method (i.e. an identical excitation 
tempera ture for all l eve l populations) would overestimate column densities only by 10 ~ 30 


percent (IHarjunpaa et ah 


20041 ) compared to the non-LTE calculation. 


The column density for those optically thin lines with one excitati on te mperature 


for all level populati ons can be calculated by the following equation ( Lee et al. 


Man gum & Shirley 


2003 


201511 


at 3/c(5 ro t 

Athin = —- T~? ex P 


kTn, 


T R dv 


(5) 


87r 3 z//i 2 J 

Here Ej = hBJ(J + 1), B is the molecular rotation constant, T ex is the excitation 
temperature, /i is the electric dipole moment, and Q rot is the partition function. In this 
calculation, the beam filling factor is assumed to be unity because, according to the SCUBA 
legacy catalogues, all our sources that have their corresponding SCUBA data have effective 
radii larger than the half of our HPBW. The partition functions are calculated up to J = 
10 (Ej =w = 233.9 K ). As described above, T ex is derived from the main isotopic line by the 
equation (4). For the sources where the main isotopic lines were not observed, we adopt the 
mean excitation temperature for the evolutionary phase that the sources belong to. 


We convert the column densities of H 13 CN and HN 13 C into those of HCN and 
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HNC us ing the fo llowing isotope ratio of carbon where -Dgc is the galactocentric radius 

~2l 


flWilson fc Root 


1994). 


12 C/ 13 C = (7.5 ± 1.9)D gc + (7.6 ± 12.9) 


( 6 ) 


To get the abundances of HCN and HNC, we also calculate the column density of molecular 
hydrogen using the (sub)millimeter continuum data. Since the majority of our HMPO 
sources were chosen from MAMBO 1.2 mm continuum survey sources, almost all of our 
HMPOs have corresponding MAMBO continuum maps. In addition, SCUBA continuum 
sources corresponding to our 4 qIRDCc, 11 aIRDCc, 18 HMPOs, and 23 UCHIIs are found. 
We tested 12 HMPO sources, which have both SCUBA and MAMBO maps and found that 
the H 2 column densities derived from different datasets are consistent with each other. 


Prior to calculating the column density with the peak flux, the dust continuum maps 
are convolved with the KVN beam (32 "). Because the dust continuum emission at 850 //m 
and 1.2 mm is optically thin, the column density of H 2 can be calculated by 


» = _ El _ 

where Tb is dust temperature, O is beam solid angle in Sr, /i is the molecul ar we ig 


(7) 


it 


ofMhe interstellar medium, and mn is the weight of one hydrogen atom fILee et al. 


20031 ). k, u is the absorption coefficient at a g iven wave length. We adopted the dust 
opacity n = 0.0102 cm 2 Q ~ 1 at 1.2 mm flKauffmann et al.l 120081 1 and n = 0.018 cm 2 g ~ 1 at 


850 /mi (iLee et al. 


2003! ). The dust temperatures are assumed to be 24 K, 35 K, 50 K and 


100 K for qIRDCc, aIRDCc, HMPOs, and UCHIIs, respectively, based on th e prev i ous 


Rathborne et al. 


studies. Dust temperatures for qIRDCc and aIRDCc are adopted from 
1 2010 ). For HMP Os, the typical Tb of ~50 K, which was estimated from the SED fitting 


bZ 


Sridharan et al 


Cesaroni et al 


( 20021) . is assumed. For UCHIIs, we adopt the gas temperature that 


(119921) have calculated using observed inversion transitions of ammonia 
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from (1,1) t° (5,5). The derived column densities, abundances, and abundance ratios are 
summarized in Tables [TIHIU 


Discussion 


4.1. Column densities and HCN/HNC abundance ratio 

Figure [4] compares the two fractional abundances of HCN and HNC. In the evolution 
from qIRDCc to UCHIIs, only the abundance of HCN increases, resulting in t 


tie increase 


of the HCN/HNC abundance ratio. This result is consistent with 


Sakai ct al 


( 2010 ), 


where no systematic variation was found in the HN 13 C abundance between the MSX 


sources (HM 
the result of 


Os) and t he MSX dark sources (IRDCs). However, this trend is opposite to 


Hirota et al 


(119981 ). where the HCN abundance is almost constant while that 


of HNC decreases as the gas temperature increases. We will discuss about this discrepancy 
in §4.2. However, we must call attention to the fact that the abundances of species highly 
depend on the adopted temperatures, which have great uncertainties. 

Table [I5l lists the mean value and standard deviation of the HCN/HNC abundance 
ratio for each evolutionary stage. The mean value gradually increases from 0.97 (± 0.10), 
2.65 (± 0.88), 4.17 (± 1.03) to 8.96 (± 3.32) in the evolutionary sequence from qIRDCc 
to UCHIIs. This indicates an increase of temperature with evolution. Figure [5] shows the 
dependence of the HCN/HNC abundance ratio on the kinetic temperature. The number of 
our sources that have the known kinetic temperature information is limited. The kinetic 
temperatures adopted for IRDCs and H 
transitions of ammonia (1,1) and (2,2) 


VIPOs have 

Deen 

Chira et al. 

2013 


Molinari et al. 


1996 ) while 


those values for UCHIIs have been calculated using the observed inversion transitions of 
ammonia from (1,1) to (5,5). A similar dependence of the HCN/HNC ratio on the kinetic 
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temperature is seen in Figure 4 of lHirota et al.l (119981 ). which claims that this dependence 
appears above 24 K. 

We also note that the mean HCN/HNC abundance ratio of aIRDCc more coincide 
with the value of HMPOs rather than that of qIRDCc within the errors. The similarity 
between aIRDCc and HMPOs is more obvious in Figure [HI which presents the cumulative 
distribution of the HCN/HNC abundance ratio for each stage. Unlike Figure [4], this 
figure represents all the objects detected both in the H 13 CN and HN 13 C lines, because the 
HCN/HNC abundance ratio only requires the column densities (or integrated intensities) 
of H 13 CN and HN 13 C. In this figure, the distribution of the HCN/HNC abundance 
ratio stretches toward a higher value as the objects evolve, but the two distributions 
for aIRDCc and HMPOs are alike. This result implies an analogous chemistry between 


the two sub-sampl es, w 


rich supports the suggestion that they are all in protostellar 


phases ({Chambers ep al. 


2009 


Sridharan et al. 


2002 


Molinari et al. 


19961 ). However, the 


mean HCN/HNC abundance ratio of HMPO is still higher than that of aIRDCc. This can 
suggest that HMPOs are relatively more evolved than aIRDCc, although both sub-samples 


are not (1 S im on et_al 


s. This suggestion is consist 

ent with the fac 

that HMPOs are 

(Ramesh & Sridharan 

1997; 

Sridharan et al. 

2002 

), but aIRDCc 


2006a) 


The separation bet ween qIRDCc and aIR DCc, however, is not obvious in the ammonia 


line surveys of IRDCs. 


Chira et al. 


(120131 ) and 


Ragan et al. 


(20111) determined the kinetic 


temperatures of the IRDC cores using the (1,1) and (2,2) ammonia inversion lines and 
reported that there is no statistically significant difference in the derived kinetic temperature 
between quiescent and active cores. According to our analysis, the two types of IRDC 
cores have a clear difference in the HCN/HNC abundance ratio, which can discriminate 
between gas temperatures, i.e., the aIRDCc have higher temperatures than the qIRDCc. 
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This discrepancy between previous work and this work may be caused by the different 
critical densities of the observed lines. The (1,1) and (2,2) ammonia lines may trace 
the cooler and more extended outer envelope than the HCN J = 1-0 line because the 
ammonia li nes have lower critical density than that of the HCN line by ~3 orders of 


magnitude flBergin fe Tafalla 


2007 


Hirota et ah 


19981 1. 


There have also been several attempts to examine the chemical and physical properties 


of massive star 


Vasvunina et al. 


brining regions using the HCN/HNC abundance ratio. For example, 


020111) studied sub-evolutionary stages of IRDC including qIRDCc and 
aIRDCc based on the 3 mm molecular line data, which includes HCN, H 13 CN and HNC 
J = 1-0 lines. In their study, kinetic temperatures (10 K ~ 30 K) that come from (1,1) 
and (2,2) ammonia inversion line surveys were adopted as T ex in calculating the molecular 
column densities. They derived the HCN/HNC abundance ratio of ~1 in IRDC cores under 
the assumption that the HNC 1-0 line is optically thin. The result is similar to our value 
for qIRDCc, but there was no si gn o f the varia tion in the HCN/HNC abundance ratio 


between qIRDCc and aIRDCc. 


Sanhueza et al. 


( 20121 ) also studied the chemical evolution 


of sub-evolutionary phases in IRDCs. They reported no significant difference in the HNC 
abundance throughout the sub-evolutionary phases of IRDCs. We also obtained similar 
abundances of HNC between the two IRDC core phases, but the HCN abundances of 
aIRDCc were larger than that of qIRDCc, making the HCN/HNC abundance ratio of the 
two IRDC cores different (see Figure QJ. 


Beyond IRDCs, 


Hoq et al. 


(120131 ) examined the HCN and HNC line data of more 


evolved objects associated with massive star formation. According to their results, 
the intensity ratio of HCN and HNC, I(HCN)/I(HNC), increased marginally from 1.07 
(prestellar cores) to 1.64 (HII regio ns/P DR s), an d the ratios were overlapped significantly 


between evolutionary stages. 


Gerner et al 


(201411 also studied the chemical evolution from 
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IRDC cores to UCHIIs. They showed that the estimated HCN/HNC abundance ratio 
ranged from only 0.3 to 0.6. However, as seen in Figure 0 and O our results are very 
different from these previous results; distinct separations of the HCN/HNC abundance 
ratio among evolutionary stages are seen in our analyses. The main reason for the different 
results between our work and the previous studies will be the optical depths of the observed 
lines. We used the optically thin lines, H 13 CN and HN 13 C, while the previous studies 
partially or fully used the H 12 CN and HN 12 C lines, which are likely optically thick. 


4.2. Abundances of HCN and HNC versus adopted temperatures 


In Figure |4] HCN abundance varies with the evolutionary stage. However, the 
abundances are highly sensitive to the adopted temperatures. The higher the adopted 
temperature is, the higher abundance is induced. This implies that the evolutionary trend 
shown in the HCN abimdance, increasing monotonically with evolution, possibly appears 
due to the adopted temperature. Excitation temperatures adopted in calculation of HCN 
column densities grow with evolution faster than those for HNC (see Table [TTJl). Thus we 
test the HCN (HNC) abundances versus the various sets of gas temperature. 


Diverse temperature information comes both from the continuum modcling^of SEDs 


and o ther molecular line obse rvations. For example, in the case of IRDCs, 


Chira et al. 


(120131 1 and 


Ragan et al. 


(20 111) determined Ty iri o 


inversion li nes. On the other hand, 


Sakai et al. 


18 K using the (1,1) and (2,2) ammonia 


line. 


Rathborne et al. 


( 20081 1 derived T rot of ~8 K using CH 3 OH 


(120101 1 derived the dust temperatures of 24 K and 35 K for qIRDC c 


and aIRD Ce. respectively, using the SED fitting method. For HMPOs, 


Sridharan et al. 


(1200211 and 


Molinari et al. 


inversion lines. 


( 199611 showed T rot of ~23 K using the (1,1) and (2,2) ammonia 


Sridharan et al 


modeling of SED, but 


(2Q0 28) de rived the dust temperatures of 50 K using the 


Beuther et al. 


(120021 ) derived a much higher temperature of 100 K in 
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their CH 3 CN observations. In the case of UCHIIs, 


Churchwell ct al. 


^30 K using the (1,1) and (2,2) ammonia inversion lines. 


(1990 


Cesaroni et al 


derived T Lin of 


(Jl992j) observed 


(4.4) and (5,5) inversion transitions of NH 3 associated with UCHII 
that the ammonia lines originate from a medium with T^ n > 50 K. 


regions and con cluded 


Churchwell ct al. 


(1199211 


suggested higher temperatures (> 100 K) in the molecular gas associated with UCHIIs 
based on the CS and CH 3 CN observation data. 


Figure |T| compares the mean abundances of HCN and HNC calculated with various 
temperature sets. The green symbols indicate the average abundances calculated in §3.2 for 
each evolutionary stage. However, for the red and blue ones, we assumed that the excitation 
temperatures are the same as the gas temperatures (i.e., the LTE condition), and adopted 
a fixed temperature for both molecules and for each evolutionary stage. The red symbols 
show the abundances derived with higher excitation temperatures, which are the same 
as the dust temperatures adopted in §3.2, while the blue ones represent the abundances 
calculated with a set of lower excitation temperatures. The adopted temperatures for the 
blue symbols are 18 K, 18 K, 23 K, 29 K for qIRDCc, aIRDCc, HMPOs, and UCHIIs, 


respectively. T 

inversion lines ( 

rese temperatures are determined from the low transition of ammonia 

Chira e 

al. 

2013; 

Ragan et al. 

2011 ; 

Sridharan et al. 

2002 ; 

Molinari et al. 

1996; 

Churchwell et al. 

1991 

!)■ 


For the higher temperatures (red symbols), both molecular abundances rapidly increase 
as objects evolve. However, for the lower temperatures (blue symbols), the HCN abundance 
slightly inc reases but the HNC abundance decreases with evolution, showing a similar 
trend to 


Hirota et ah 


(119981) . This indicates that the abundance variation shown in our 


calculations (Figure [4]) could be significantly affected by the adopted temperatures. 


However, even though we adopt the different sets of excitation temperature, the 
HCN/HNC abundance ratio still increases from 0.83(± 0.06), 2.02(± 0.6) 2.65(± 0.5) to 
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4.9(± 1.4) with evolution (blue and red symbols in Figure [71). This show that the increase of 
the HCN/HNC abundance ratio is originated not from the adopted excitation temperatures 
but from the evolution of kinetic environment. 


4.3. Chemical Modeling 


The variation of the HCN/HNC abundance ratio with temperature has been a debated 
issue. The HCN/HNC ab undance ratio has been calcul ated in various sou rces, from dark 


clouds (IHirota et al. 


19981 1 . active star-formi ng region (Schilke et ah 


galaxies and associated large-scale outflows ( Aalto et al. 


2012 ). 


1992b to external 


Schilke et al. 


(119921 ) 


suggested that the neutral-neutral reaction described in Equation (2) mainly governs the 
transition of the HCN/HNC abundance ratio, and set the activation energy of the reaction 


to be 200 K, which best reprodu ces t he HCN 
determined. On the other hand, 


Talbietal. 


/ HNC abundance ratio observationally 


199611 quantum-chemically analyzed the 


neutral-neutral reaction and suggested that the reaction should have an activation energy 
as high as 2000 K, resulting in a rate coefficient of ~ 10 “ 15 at a temperature lower 
than 100 K. However, this rate cannot reproduce the observed results. The discrepancy 
between the rate coefficients empirically determined (h ereafter k^rh^) and t heoretically 


determined (hereafter h Ta ibi) has remained unresolved (Graninger et al. 


201411. 


We tested both rate coefficients for our chemical calculations and compared them 
with our observations. We modeled the HCN/H NC ab undance ratio evolution using the 


astrochemical cod e, ALCHEMIC (ISemenov et al. 


reaction network ( McElrov et al 


2010 1 with the UMIST gas chemical 


201311 . The only surface chemistry considered in 


our calculation is t 


Cazaux fe Tielens 


re H -2 formation, and its rate equation has been adopted from 


2004, 


20101 1. The chemical evolution has been calculated adopting 


three continuous evolutionary stages of physical conditions: from IRDCs, through HMPOs, 
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to UCHIIs. The physical conditions at a given evolutionary stage are assumed constant. 
The initial chemical abundances for the HMPO stage were adopted from the final chemical 
abundances for the IRDC stage, and the final chemical abundances for the stages of HMPOs 
were taken as the initial abundances for the UCHII stage. The initial abu ndance s for the 


Gerner et al. 


( 120 1 41 )) 


IRDC stage were adopted from the set with low metals (Table 7 of 
However, the atomic abundance of N is decreased from 2.47 x 10 -5 to 1.00 x 10 -5 to achieve 
a better fit to the calculated abundances. The physical condition and chemical timescale 
adopted for each phase are listed in Table [161 Since the chemical evolution is associated 
with kinetic temperature, in our chemical calculations, we assume that the gas temperatures 
are the same as the dust temperatures adopted in §3.1. In this dense environment, gas and 
dust are thermally well coupled. Actually we tested a variety of gas temperature in our 
chemical calculation. According to our test, lower gas temperatures adopted for the blue 
symbols in Figure [7] cannot reproduce the observed HCN/HNC abundance ratios. 

The model mainly solves the following equation for 4604 gas-phase reactions. 

d n(i) 


d£ 


^ ^ Tli ^ ^ 

i^l 


( 8 ) 


l.m 


where n(i) is the gas-phase number density of the i-th species (cm 3 ), and ki m and ki 


are rates of the gas-phase reactions. For the first-order kin etics and for_t 
kinetics, the units are s -1 and cm 3 s _1 , respectively ( Semenov et al. 


le second-order 


20101). The cosmic 


ray (CR) ionization rate is assumed to be 5.0 x 10 -17 s” 1 . For simplicity, the effect by UV 
photons was not included, i.e., the model gas is assumed to be well shielded from either 
external or internal UV radiation. This assumption could be reasonable because HCN 
traces very dense gas; the dense inner region of prestellar cores will be shielded from the 
external interstellar radiation held, which is the only UV radiation source in this stage. 
On the other hand, in the protostellar stage, the UV photons emitted from the central 
protostars will be absorbed by the very dense infalling material, which prevents the high 
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energy photons from traveling farther out. In addition, the UV photons from the early type 
central stars in UCHIIs are absorbed by the inne r warm dusty en velope, which converts the 


entire stellar luminosity to FIR radiation (jChurchwell. Ed. 


2002 ). 


Figure 0 compares the observed HCN/HNC abund ance ratio with the calculated 


values. Here the chemical time scales are taken from 


Gerner et al. 


(120141 1. Different line 


types indicate different activation energies used in the chemical model. For the observed 
HCN/HNC abundance ratio of IRDCs, we adopt the mean ratio of all qIRDCc and aIRDCc 
since the timescale for each of those cores is not well known. The calculated HCN/HNC 
abundance ratio is controlled by the activation energy for the neutral-neutral reaction, 
and that fcs c hiike provides results much more consistent with our observational ratios 
compared to A/paibi- However, the observed ratios are better fitted by the model with an 
activation energy of 150 K, which results in a rate coefficient higher than &s c hiike- In this 
chemical calculations, we did not try to fit the exact abundances of HCN and HNC since 
the abundances calculated from observations are very sensitive to the adopted excitation 
temperature as mentioned in the previous section. 


The model we adopt here is too simple to provide a precise analysis; only three discrete 
physical conditions have been adopted to describe the evolution of massive star formation, 
and the timescale for each evolutionary stage has large uncertainty. Nevertheless, our 
test with chemical models clearly show that (1) the HCN/HNC abundance ratio can be 
a evolutionary tracer of massive star formation and (2) to fit the ratios derived from 
observations, the activation energy for the neutral-neutral reaction should be much smaller 
than the theoretical value. 
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5. Summary 

To examine evolutionary effects in the chemical and physical conditions of massive 
star-forming regions, we carried out H 13 CN and HN 13 C line observations and analyzed 
the HCN/HNC abundance ratio toward IRDC cores (qIRDCc and aIRDCc), HMPOs, and 
UCHIIs, which are known as early phases of massive star formation. We detected both 
H 13 CN and HN 13 C lines toward 8 qIRDCc, 16 aIRDCc, 23 HMPOs, and 31 UCHIIs. 

1. We found a statistically increasing tendency of the abundance ratio with the 
evolution of the objects, from 0.97 (qIRDCc) to 8.96 (UCHIIs). This follows the known 
evolutionary scheme for massive star formation well and indicates that the HCN/HNC 
abundance ratio can be used to trace the evolutionary stages of massive star formation, in 
which temperature increases. 

2. aIRDCc and HMPOs show similar mean values and distributions of the HCN/HNC 
abundance ratio. This implies that the two types of sources have similar chemical conditions 
and supports the suggestion that both types of sources may be in the protostellar phase. 

3. Although aIRDCc and HMPO show similar abundance ratios, the mean value of 
HMPOs is still higher than that of aIRDCc. This may indicate that both sub-samples are 
in the same protostellar phase, but HMPOs are relatively more evolved than aIRDCc. This 
suggestion is consistent with the fact that HMPOs are associated with MSX sources, while 
aIRDCc are not. 

4. Our observed HCN/HNC abundance ratio cannot be explained by the theoretically 
determined rate coefficients of the neutral-neutral reaction. This discrepancy between 
observational and theoretical works remains to be resolved. 
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Fig. 1.— The representative observed spectra for four sub-samples. The vertical dotted line 
indicates the peak velocity of the HN 13 C J=l-0 line for each source (Tables 6-9). 
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Fig. 2.— The averaged H 13 CN and HN 13 C spectra for each sub-sample. The x-axis is 
relative velocity with respect to the HN 13 C peak velocity of each source, while the y-axis is 
beam corrected antenna temperature. 
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Fig. 3.— Comparisons of the line widths between H 13 CN and HN 13 C. The dotted lines 
indicate 1-cr for the linear regression fits. The Pearson correlation coefficient (R) and p- 
value are given in the upper left corner of each box. 
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O 



X 



X[HNC] 


Fig. 4.— Plot of the abundance of HCN against that of HNC. Only the sources which have 
their corresponding continuum map are plotted. 1RDC cores, HMPOs, and UCHIIs are 
indicated by diamonds, times, and crosses, respectively. Here IRDC cores are divided into 
quiescent ones (filled diamonds) and active ones (open diamonds). Observational uncertainty 
is smaller than the symbol size. 
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Fig. 5.— Plot of the HCN/HNC abundance ratio versus the kinetic temperature for our 


sourc es that have t heir k inetic tempe rature informat i on. T 
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( 2013 ) for IRDCc, 


Molinari et al 


le kinetic tempera 


lure is adopted 


( 19961 ) for HMPO, and 


Cesaroni et ah 


























32 


Cumulative Distribution 



Fig. 6.— Cumulative distribution of the HCN/HNC abundance ratio. qIRDCc and UCHIIs 
show clearly different distributions from the other sub-samples. On the contrary, aIRDCc 
and HMPOs show similar distributions, implying the similar chemical conditions between 
the two sub-samples. 
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X[HNC] 


Fig. 7.— The mean abundance of HCN against that of HNC derived from different sets of 
temperature. All symbols are the same as in Figure^ The green symbols indicate the average 
abundances calculated with different excitation temperatures for H 13 CN and HN 13 C (see 
§3.2) for each evolutionary stage. However, for the red (the highest temperature set) and 
blue (the lowest temperature set) ones, we assumed the same excitation temperatures as the 
dust temperatures, which are fixed for each evolutionary stage. See the text for details. 
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Fig. 8.— A comparison between the HCN/HNC abundance ratio derived from the ob¬ 
servation and the ratios calculated from the chemical model. Dashed, dotted, an d solid 


lines indicate the cai 


200 K ( Schilke et al. 


culated ratios with activation barriers of 2000 K flTalbi et ah 
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Table 1. Observed Lines 


Molecule 

Transition 

v [MHz] 

MD] 

Ca 

°ul 

HCN 

J=l-0, F=2-l 

88630.42 

2.99 b 

3 


J=l-0, F=l- 1 

88631.85 


5 


J=l-0, F=0-1 

88633.94 


1 

h 13 cn 

J=l-0, F=2-l 

86338.77 

2.99 b 

3 


J=l-0, F=l-1 

86340.18 


5 


J=l-0, F=0-1 

86342.27 


1 

HNC 

J=l-0 

90663.57 

3.05 c 


hn 13 c 

J=l-0 

87090.85 

3.05 c 



Note. _ — a In tr insic line strength. 


Bhattacharva & Gordv 

( 

1960 

)• c 

Blackman et al. 


1976|). 
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Table 2. 

. qIRDCc 

source 

Information 



Source Name 

RA 

J(2000.0) 

Dec 

J(2000.0) 

l 

[°] 

b 

[°] 

Dkin 

[kpc] 

Dgc 

[kpc] 

R 

["] 

Tgas 

[K\ 

G018.82_MM6 

18:26:18.4 

-12:41:15 

18.80 

-0.29 

4.8 

5.7 



G028.37JVIM9 

18:42:46.7 

-04:04:08 

28.32 

0.07 

5.0 

5.3 

57.1 


G028.37JVIM17 

18:43:00.0 

-04:01:34 

28.39 

0.04 

5.0 

5.3 

40.7 

9.3 

G028.53JVIM7 

18:44:23.7 

-04:02:09 

28.54 

-0.28 

5.7 

5.1 


11.7 

G031.97_MM9 

18:49:31.6 

-00:46:30 

32.02 

0.07 

6.9 

5.1 



G033.69_MM8 

18:52:53.9 

+00:41:16 

33.71 

-0.01 

7.1 

5.0 

27.9 


G034.77JV1M4 

18:56:48.9 

+01:23:34 

34.78 

-0.56 

2.9 

6.7 


16.3 

G035.39JV1M5 

18:57:08.8 

+02:08:09 

35.48 

-0.30 

2.9 

6.7 


19.7 


Note. — Kinetic distances are quoted from 


Simon et al 


( 2006b '). The source 


(Di Francesco et al. 

2008 

). The Tgas are quoted from Tki n in 

Chira et al. 

(2013) 
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Tabic 3. aIRDCc source Information 


Source Name 

RA 

J(2000.0) 

Dec 

J(2000.0) 

l 

[°] 

b 

[°] 

-Dkin 

[kpc] 

Dgc 

[kpc] 

Ii 

H 

Tgas 

[K] 

G019.27_MM2 

18:25:52.6 

-12:04:48 

19.29 

0.08 

2.4 

6.6 

41.0 

9.2 

G022.35_MM1 

18:30:24.4 

-09:10:34 

22.38 

0.45 

4.3 

5.7 

58.0 

15.8 

G023.60_MM1 

18:34:11.6 

-08:19:06 

23.57 

0.01 

3.9 

4.3 



G024.60_MM1 

18:35:40.2 

-07:18:37 

24.63 

0.15 

3.7 

5.9 

49.3 

20.2 

G028.37_MM4 

18:42:50.7 

-04:03:15 

28.34 

0.06 

5.0 

5.0 

51.0 

13.0 

G028.37_MM6 

18:42:49.0 

-04:02:23 

28.35 

0.07 

5.0 

5.3 

40.9 


G030.97_MM1 

18:48:21.6 

-01:48:27 

30.97 

-0.14 

5.1 

5.4 

54.8 


G031.97JVIM1 

18:49:36.3 

-00:45:45 

32.04 

0.06 

6.9 

5.1 

53.3 


G033.69_MM4 

18:52:56.4 

+00:43:08 

33.74 

-0.01 

7.1 

5.0 

44.4 

13.1 

G033.69_MM5 

18:52:47.8 

+00:36:47 

33.63 

-0.02 

7.1 

5.0 

26.6 


G034.43_MM1 

18:53:18.0 

+01:25:24 

34.41 

0.23 

3.7 

6.2 

39.6 


G034.43JV1M3 

18:53:20.4 

+01:28:23 

34.46 

0.25 

3.7 

6.2 

43.4 


G034.43_MM4 

18:53:19.0 

+01:24:08 

34.39 

0.22 

3.7 

6.2 



G035.39_MM7 

18:57:08.1 

+02:10:50 

35.52 

-0.27 

2.9 

6.7 



G048.65_MM1 

19:21:49.7 

+13:49:30 

48.67 

-0.30 

2.5 

7.4 

26.6 


G053.25_MM6 

19:29:31.5 

+17:59:50 

53.22 

0.05 

1.9 

7.8 




Note. 


Kinetic distances are quoted from 


Simon et al 


(|2006bh .Tlie source 


sizes (R) are quoted fro m effective radius in the JCMT SCUBA legacy c atalogue 


(Di Francesco et al. 

2008 

). The T gas are quoted from Tk; n in 

Chira et al. 

(2013) 















Table 4. HMPO source Information 


Source Name 

RA 

J(2000.0) 

Dec 

J(2000.0) 

l 

[°] 

b 

[1 

-£^kin,far 

[kpc] 

-£^kin,near 

[kpc] 

Dgc 

[kpc] 

R 

["] 

Tgas 

\K] 

IRAS05358+3543 

05:39:10.4 

+35:45:19 

173.48 

2.43 

1.8 


21.2 

47.3 

18.0 b 

IRAS05373+2349 

05:40:24.4 

+23:50:54 

183.72 

-3.66 

1.2 


9.8 

42.5 

21.2“ 

IRAS18024-2119 

18:05:25.4 

-21:19:41 

8.83 

-0.03 

0.1 


8.4 



IRAS18089-1732 

18:11:51.3 

-17:31:28 

12.89 

0.49 

13.0 

3.6 

5.5 

44.9 

19.2“ 

IRAS18102-1800 

18:13:12.2 

-17:59:35 

12.63 

-0.02 

14.0 

2.6 

6.3 


15.0 b 

IRAS18144-1723 

18:17:24.4 

-17:22:13 

13.66 

-0.60 

4.3 


4.9 

58.3 

23.6“ 

IRAS18151-1208 

18:17:57.1 

-12:07:22 

18.34 

1.77 

3.0 


6.2 


23.7“ 

IRAS18182-1433 

18:21:07.9 

-14:31:53 

16.58 

-0.05 

11.8 

4.5 

4.8 

38.0 

20.0 b 

IRAS18223-1243 

18:25:10.9 

-12:42:17 

18.66 

-0.06 

12.4 

3.7 

5.6 


18.0 b 

IRAS18264-1152 

18:29:14.3 

-11:50:26 

19.88 

-0.53 

12.5 

3.5 

5.8 

24.9 

18.0 b 

IRAS18308-0841 

18:33:31.9 

-08:39:17 

23.20 

0.00 

10.7 

4.9 

4.9 

30.5 

18.0 b 

IRAS18345-0641 

18:37:16.8 

-06:38:32 

25.41 

0.10 

9.5 


4.6 

42.4 

16.0 b 

IRAS18440-0148 

18:46:36.3 

-01:45:23 

30.82 

0.27 

8.3 


5.0 

30.5 

23.0 b 

IRAS18470-0044 

18:49:36.7 

-00:41:05 

32.11 

0.09 

8.2 


5.1 


20.0 b 

IRAS18511+0146 

18:53:38.1 

+01:50:27 

34.82 

0.35 

3.9 


6.2 

60.5 

95.5“ 

IRAS18530+0215 

18:55:34.2 

+02:19:08 

35.47 

0.14 

8.7 

5.1 

5.7 


16.0 b 

IRAS19220+1432 

19:24:19.7 

+14:38:03 

49.67 

-0.46 

5.5 


6.4 

21.4 

23.0 b 

IRAS19410+2336 

19:43:11.4 

+23:44:06 

59.78 

0.06 

6.4 

2.1 

7.9 

51.2 

18.0 b 

IRAS19411+2306 

19:43:18.1 

+23:13:59 

59.36 

-0.21 

5.8 

2.9 

7.7 

37.1 

14.0 b 

IRAS20126+4104 

20:14:26.0 

+41:13:32 

78.12 

3.63 

1.7 


8.8 

47.5 

28.1“ 

IRAS20343+4129 

20:36:07.1 

+41:40:01 

80.83 

0.57 

1.4 


8.3 

61.5 

18.0 b 

IRAS22198+6336 

22:21:27.6 

+63:51:42 

107.30 

5.64 

1.3 


9.0 


21.5“ 

IRAS23033+5951 

23:05:25.7 

+60:08:08 

110.09 

-0.07 

3.5 


11.0 

47.9 

20.0 b 
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Note. — Kinetic distances are quoted from 


Beuthcr et al 


( 2002 ) and from 


Molinari et al. 


( 19961) . If the 


distance ambiguity is resolved, only far and no near dista nce is noted. The source s izes ( R ) are quoted from 


from Tkm i n 


Molinari et al. 

(1996 

) and the The T gas b are quoted from T rot in 

Beuther et al. 

(2002) 
























Table 5. UCHII source Information 


Source Name 

RA 

J(2000.0) 

Dec 

J(2000.0) 

l 

[°] 

b 

[°] 

-Dkin 

[kpc] 

Dgc 

[kpc] 

R 

["] 

Tgas 

[K\ 

IRAS02232+6138 

02:27:01.0 

+61:52:14 

133.94 

1.06 

3.0 

11.4 



IRAS02575+6017 

03:01:32.3 

+60:29:12 

138.30 

1.56 

3.8 

11.0 

62.0 


IRAS03035+5819 

03:07:25.6 

+58:30:52 

139.91 

0.20 

4.2 

11.2 

48.0 


IRAS06053-0622 

06:07:46.6 

-06:22:59 

213.70 

-12.60 

10.8 

9.2 



IRAS06084-0611 

06:10:51.0 

-06:11:54 

213.88 

-11.84 

1.0 

9.3 



IRAS06056+2131 

06:08:41.0 

+21:31:01 

189.03 

0.78 

0.8 

8.6 

63.4 


IRAS06058+2138 

06:08:54.1 

+21:38:25 

188.95 

0.89 

2.2 

9.2 



IRAS06061+2151 

06:09:07.8 

+21:50:39 

188.80 

1.03 

4.1 

8.5 

59.4 


IRAS06099+1800 

06:12:53.3 

+17:59:22 

192.60 

-0.05 

2.5 

10.4 

45.8 


IRAS17574-2403 

18:00:30.4 

-24:04:00 

5.89 

-0.39 

2.0 

6.3 

33.5 

200.0 

IRAS17599-2148 

18:03:00.4 

-21:48:05 

8.14 

0.23 

4.2 

5.6 



IRAS18032-2137 

18:06:19.0 

-21:37:32 

8.67 

-0.36 

4.8 

4.4 

44.6 


IRAS18032-2032 

18:06:13.5 

-20:31:47 

8.67 

-0.36 

5.7 

7.7 


129.0 

IRAS18075-1956 

18:10:23.5 

-19:56:15 

10.61 

-0.37 

4.8 

9.2 



IRAS18162-2048 

18:19:11.9 

-20:47:34 

10.84 

-2.59 

1.9 

6.1 



IRAS18100-1854 

18:14:01.1 

-18:53:24 

11.94 

-0.62 

5.2 

5.0 



IRAS18174-1612 

18:20:24.8 

-16:11:35 

15.03 

-0.68 

2.1 

6.6 



IRAS18317-0757 

18:34:24.9 

-07:54:48 

23.95 

0.15 

6.0 

4.9 

41.2 


IRAS18316-0602 

18:34:19.8 

-05:59:44 

25.65 

1.05 

3.0 

6.3 



IRAS18403-0417 

18:42:58.2 

-04:14:00 

28.20 

-0.05 

9.1 

4.8 

47.5 


IRAS18479-0005 

18:50:30.9 

-00:01:59 

32.80 

0.19 

12.9 

7.8 



IRAS18507+0110 

18:53:18.5 

+01:14:58 

34.26 

0.15 

4.0 

6.2 

42.5 


IRAS19095+0930 

19:11:53.3 

+09:35:46 

43.79 

-0.13 

9.0 

6.9 

56.2 


IRAS19110+1045 

19:13:22.1 

+10:50:53 

45.07 

0.13 

6.0 

6.6 

29.7 









Table 5—Continued 


Source Name 

RA 

J(2000.0) 

Dec 

J(2000.0) 

l 

[°] 

b 

[°] 

-Dkin 

[kpc] 

Dgc 

[kpc] 

R 

["] 

Tgas 

[K] 

IRAS20081+3122 

20:10:09.1 

+31:31:34 

69.54 

-0.98 

3.0 

8.2 



IRAS20255+3712 

20:27:26.6 

+37:22:48 

76.38 

-0.62 

1.0 

8.7 

42.3 


IRAS20178+4046 

20:19:39.3 

+40:56:30 

78.44 

2.66 

3.3 

8.6 



IRAS20350+4126 

20:36:52.6 

+41:36:32 

80.87 

0.42 

2.1 

8.7 

63.4 


IRAS22543+6145 

22:56:19.1 

+62:01:57 

109.87 

2.12 

0.7 

9.0 



IRAS23138+5945 

23:16:04.8 

+60:02:00 

111.28 

-0.66 

2.5 

10.5 

59.4 


IRAS23133+6050 

23:15:31.5 

+61:07:09 

111.61 

0.37 

5.2 

11.2 

70.5 



Note. — Kinetic distances are quoted from 


Thompson et al 


The source sizes [R 

are c 

(Di Francesco et al. 

2008) 


200811 . The T gas are quoted from Tki n in 


(12006T ) and references therein. 

catalogue 


Cesaroni et al 


( 199211 . 
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Table 6. qIRDCc Observed Line Parameters 


Source Name 

Molecule 

Transition 

TmB 

[K] 

Trms 

[K] 

Av 

[kms -1 ] 

/T MB dv 
[K kms -1 ] 

■ULSR 

[kms -1 ] 

T thin 

Tex 

[K] 

G018.82.MM6 

h 13 cn 

F = 0-1 

0.039 

0.034 

2.407 ( 1.710 ) 

0.099 ( 

; 0.054 ) 

65.7 





F = 2-1 

0.132 

0.034 

1.302 ( 0.410 ) 

0.183 ( 

; 0.067 ) 






F = 1-1 

0.085 

0.034 

6.946 ( 1.220 ) 

0.626 ( 

; o.ii8) 





hn 13 c 


0.544 

0.057 

1.751 ( 0.139 ) 

1.014 ( 

; 0.073 ) 




G028.37.MM9 

h 13 cn 

F = 0-1 

0.047 

0.047 

3.222 ( 1.610 ) 

0.162 ( 

; 0.078 ) 

80.1 

0.13 

12.55 



F = 2-1 

0.190 

0.047 

3.267 ( 0.487 ) 

0.662 ( 

; 0.088 ) 






F = 1-1 

0.141 

0.047 

3.457 ( 0.691 ) 

0.518 ( 

; 0.091) 





hn 13 c 


0.625 

0.060 

2.548 ( 0.157 ) 

1.695 ( 

; 0.093 ) 


0.14 

13.39 

G028.37.MM17 

h 13 cn 

F = 0-1 

0.044 

0.051 

17.690 ( 4.720 ) 

0.820 ( 

; 0.281 ) 

79.1 





F = 2-1 

0.210 

0.051 

1.839 ( 0.549 ) 

0.410 ( 

: o.in) 






F = 1-1 

0.090 

0.051 

2.096 ( 1.020 ) 

0.200 ( 

; 0.105) 





hn 13 c 


0.509 

0.057 

2.631 ( 0.207 ) 

1.424 ( 

; 0.093 ) 




G028.53.MM7 

h 13 cn 

F = 0-1 

0.082 

0.030 

2.054 ( 1.060 ) 

0.179 ( 

; 0.057 ) 

88.5 





F = 2-1 

0.146 

0.030 

1.613 ( 0.269 ) 

0.250 ( 

; 0.038 ) 






F = 1-1 

0.101 

0.030 

1.747 ( 0.663 ) 

0.188 ( 

; 0.047 ) 





hn 13 c 


0.588 

0.063 

1.637 ( 0.172 ) 

1.024 ( 

; 0.084 ) 




G031.97.MM9 

h 13 cn 

F = 0-1 






96.8 

0.13 

28.79 



F = 2-1 











F = 1 - 1 










hn 13 c 


0.243 

0.045 

5.556 ( 0.465 ) 

1.436 ( 0.105 ) 


0.09 

18.07 

G033.69.MM8 

h 13 cn 

F = 0 - 1 






103.3 





F = 2 - 1 











F = 1 - 1 










hn 13 c 


0.127 

0.025 

7.385 ( 0.517 ) 

0.999 ( 0.065 ) 




G034.77.MM4 

h 13 cn 

F = 0 - 1 

0.031 

0.040 

8.712 ( 4.350 ) 

0.288 ( 0.120 ) 

41.7 





F = 2 - 1 

0.196 

0.040 

1.513 ( 0.279 ) 

0.315 ( 0.047 ) 






F = 1 - 1 

0.110 

0.040 

1.679 ( 0.361 ) 

0.196 ( 0.047 ) 





hn 13 c 


0.531 

0.042 

1.840 ( 0.115 ) 

1.041 ( 0.057 ) 




G035.39.MM5 

h 13 cn 

F = 0 - 1 

0.071 

0.047 

0.434 ( 1.070 ) 

0.033 ( 0.030 ) 

45.5 

0.07 

15.04 



F = 2 - 1 

0.201 

0.047 

2.420 ( 0.420 ) 

0.517 ( 0.074 ) 






F = 1 - 1 

0.149 

0.047 

2.015 ( 0.536 ) 

0.320 ( 0.068 ) 





hn 13 c 


0.529 

0.058 

1.835 ( 0.216 ) 

1.033 ( 0.088 ) 


0.11 

13.85 
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Tabic 7. AIRDC Observed Line Parameters 


Source Name 


G019.27.MM2 


G022.35.MM1 


G023.60.MM1 


G024.60.MM1 


G028.37.MM4 


G028.37.MM6 


G030.97.MM1 


G031.97.MM1 


G033.69.MM4 


G033.69.MM5 


Molecule 

Transition 

t mb 

[K] 

^rms 

[K] 

Av 

[kms -1 ] 

J r MB dv 
[K kms -1 ] 

^LSR 

[kms -1 ] 

^thin 

Tex 

[K] 

h 13 cn 

F = 0-1 

0.104 

0.052 

1.881 ( 0.746 ) 

0.209 ( 0.072 ) 

26.7 

0.16 

10.22 


F = 2-1 

0.263 

0.052 

2.388 ( 0.393 ) 

0.668 ( 0.091 ) 





F = 1-1 

0.231 

0.052 

2.315 ( 0.722 ) 

0.570 ( 0.115 ) 




hn 13 c 


0.497 

0.069 

1.578 ( 0.214 ) 

0.834 ( 0.090 ) 


0.18 

7.00 

h 13 cn 

F = 0-1 

0.145 

0.063 

1.488 ( 0.439 ) 

0.229 ( 0.071 ) 

52.6 

0.16 

19.70 


F = 2-1 

0.311 

0.063 

2.988 ( 0.427 ) 

0.990 ( 0.115 ) 





F = 1-1 

0.203 

0.063 

3.492 ( 0.789 ) 

0.756 ( 0.130 ) 




hn 13 c 


0.391 

0.049 

2.191 ( 0.248 ) 

0.913 ( 0.078 ) 




h 13 cn 

F = 0-1 





107.3 

0.30 

36.93 


F = 2-1 









F = 1-1 








hn 13 c 


0.326 

0.066 

8.984 ( 2.020 ) 

3.121 ( 0.389 ) 


0.30 

16.35 

h 13 cn 

F = 0-1 

0.050 

0.061 

7.623 ( 2.870 ) 

0.410 ( 0.165 ) 

53.1 

0.16 

16.77 


F = 2-1 

0.258 

0.061 

5.834 ( 0.907 ) 

1.604 ( 0.197 ) 





F = 1-1 

0.149 

0.061 

2.172 ( 0.737 ) 

0.346 ( 0.130 ) 




hn 13 c 


0.379 

0.055 

1.890 ( 0.199 ) 

0.763 ( 0.073 ) 


0.06 

11.80 

h 13 cn 

F = 0 - 1 

0.060 

0.053 

8.614 ( 3.950 ) 

0.549 ( 0.207 ) 

79.2 

0.16 

20.76 


F = 2 - 1 

0.420 

0.053 

3.062 ( 0.348 ) 

1.368 ( 0.180 ) 





F = 1 - 1 

0.301 

0.053 

3.676 ( 0.574 ) 

1.179 ( 0.144 ) 




hn 13 c 


0.757 

0.073 

3.009 ( 0.190 ) 

2.426 ( 0.128 ) 


0.15 

19.23 

h 13 cn 

F = 0 - 1 

0.166 

0.050 

2.406 ( 0.609 ) 

0.426 ( 0.083 ) 

80.3 

0.30 

11.95 


F = 2 - 1 

0.358 

0.050 

3.356 ( 0.321 ) 

1.281 ( 0.149 ) 





F = 1-1 

0.181 

0.050 

5.685 ( 1.090 ) 

1.094 ( 0.177 ) 




hn 13 c 


0.839 

0.053 

2.407 ( 0.109 ) 

2.150 ( 0.082 ) 


0.27 

10.20 

h 13 cn 

F = 0 - 1 

0.131 

0.042 

3.025 ( 0.655 ) 

0.421 ( 0.074 ) 

78.0 

0.09 

25.88 


F = 2-1 

0.411 

0.042 

1.930 ( 0.139 ) 

0.845 ( 0.057 ) 





F = 1 - 1 

0.193 

0.042 

3.311 ( 0.777 ) 

0.681 ( 0.101 ) 




hn 13 c 


0.462 

0.047 

2.549 ( 0.191 ) 

1.254 ( 0.077 ) 


0.09 

17.25 

h 13 cn 

F = 0 - 1 

0.270 

0.047 

4.508 ( 0.793 ) 

1.295 ( 0.183 ) 

95.4 

0.17 

54.37 


F = 2 - 1 

0.880 

0.047 

3.767 ( 0.269 ) 

3.527 ( 0.254 ) 





F = 1 - 1 

0.659 

0.047 

4.574 ( 0.361 ) 

3.211 ( 0.205 ) 




hn 13 c 


0.657 

0.043 

3.548 ( 0.138 ) 

2.482 ( 0.081 ) 


0.11 

26.90 

h 13 cn 

F = 0 - 1 

0.088 

0.054 

1.460 ( 0.724 ) 

0.137 ( 0.057 ) 

106.1 

0.13 

15.87 


F = 2 - 1 

0.223 

0.054 

3.673 ( 0.668 ) 

0.871 ( 0.126 ) 





F = 1 - 1 

0.194 

0.054 

3.233 ( 0.568 ) 

0.667 ( 0.101 ) 




hn 13 c 


0.501 

0.043 

3.342 ( 0.180 ) 

1.782 ( 0.080 ) 


0.16 

13.28 

h 13 cn 

F = 0 - 1 

0.062 

0.049 

0.664 ( 0.555 ) 

0.044 ( 0.037 ) 

105.3 

0.09 

11.31 


F = 2 - 1 

0.215 

0.049 

1.892 ( 0.372 ) 

0.433 ( 0.069 ) 





F = 1 - 1 

0.132 

0.049 

1.387 ( 0.752 ) 

0.194 ( 0.071 ) 




hn 13 c 


0.417 

0.041 

1.865 ( 0.164 ) 

0.827 ( 0.058 ) 


0.08 

12.47 

h 13 cn 

F = 0 - 1 

0.399 

0.074 

15.740 ( 1.120 ) 

6.679 ( 0.513 ) 

57.8 

0.17 

41.36 


F = 2 - 1 

0.774 

0.074 

2.747 ( 0.184 ) 

2.264 ( 0.245 ) 





G034.43.MM1 
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Table 7—Continued 


Source Name 

Molecule 

Transition 

TmB 

[K] 

^rms 

[K] 

Av 

[kms -1 ] 

f ^MB^v 
[K kms -1 ] 

^LSR 

[kms -1 ] 

T thin 

Tex 

[K] 



F = 1-1 

0.577 

0.074 

2.657 ( 0.314 ) 

1.633 ( 0.246 ) 





hn 13 c 


0.749 

0.049 

2.863 ( 0.113 ) 

2.282 ( 0.080 ) 


0.07 

31.26 

G034.43.MM3 

h 13 cn 

F = 0-1 

0.140 

0.044 

13.080 ( 3.870 ) 

1.956 ( 0.654 ) 

59.4 

0.13 

35.33 



F = 2-1 

0.400 

0.044 

2.205 ( 0.348 ) 

0.938 ( 0.218 ) 






F = 1-1 

0.216 

0.044 

4.901 ( 0.809 ) 

1.129 ( 0.351 ) 





hn 13 c 


0.790 

0.051 

2.174 ( 0.111 ) 

1.827 ( 0.076 ) 


0.10 

23.11 

G034.43.MM4 

h 13 cn 

F = 0-1 

0.199 

0.047 

4.897 ( 0.776 ) 

1.036 ( 0.128 ) 

57.8 

0.27 

25.59 



F = 2-1 

0.736 

0.047 

3.488 ( 0.223 ) 

2.733 ( 0.148 ) 






F = 1 - 1 

0.481 

0.047 

3.737 ( 0.332 ) 

1.915 ( 0.136 ) 





hn 13 c 


0.952 

0.050 

2.525 ( 0.096 ) 

2.559 ( 0.080 ) 


0.17 

18.53 

G035.39.MM7 

h 13 cn 

F = 0 - 1 

0.082 

0.033 

0.434 ( 0.845 ) 

0.038 ( 0.022 ) 

45.8 

0.09 

10.85 



F = 2-1 

0.237 

0.033 

1.864 ( 0.273 ) 

0.470 ( 0.050 ) 






F = 1 - 1 

0.122 

0.033 

2.711 ( 0.568 ) 

0.353 ( 0.057 ) 





hn 13 c 


0.671 

0.065 

1.994 ( 0.141 ) 

1.424 ( 0.090 ) 


0.16 

11.07 

G048.65.MM1 

h 13 cn 

F = 0 - 1 

0.091 

0.048 

0.434 ( 0.586 ) 

0.042 ( 0.032 ) 

34.0 

0.09 

8.73 



F = 2 - 1 

0.275 

0.048 

0.910 ( 0.185 ) 

0.266 ( 0.044 ) 






F = 1 - 1 

0.177 

0.048 

1.747 ( 0.358 ) 

0.328 ( 0.063 ) 





hn 13 c 


0.237 

0.048 

1.658 ( 0.260 ) 

0.418 ( 0.059 ) 


0.05 

8.57 

G053.25.MM6 

h 13 cn 

F = 0 - 1 

0.182 

0.041 

0.480 ( 6.180 ) 

0.093 ( 0.035 ) 

23.7 

0.04 

12.51 



F = 2 - 1 

0.194 

0.041 

1.263 ( 0.189 ) 

0.260 ( 0.041 ) 






F = 1 - 1 

0.116 

0.041 

1.063 ( 0.318 ) 

0.131 ( 0.039 ) 





hn 13 c 


0.328 

0.059 

1.445 ( 0.193 ) 

0.504 ( 0.067 ) 


0.02 

9.56 
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Tabic 8. HMPO Observed Line Parameters 


Source Name 


IRAS05358+3543 


IRAS05373+2349 


IRAS 18024-2119 


IRAS 18089-1732 


IRAS 18102-1800 


IRAS 18144-1723 


IRAS 18151-1208 


IRAS 18182-1433 


IRAS 18223-1243 


IRAS 18264-1152 


Molecule 

Transition 

t mb 

[K] 

Irms 

[K] 

Av 

[kms -1 ] 

I r MB dv 
[K kms -1 ] 

^LSR 

[kms -1 ] 

T thin 

Tex 

[K] 

h 13 cn 

F = 0-1 

0.075 

0.044 

6.992 ( 2.300 ) 

0.557 ( 

; 0.166 ) 

-17.4 

0.16 

23.50 


F = 2-1 

0.444 

0.044 

2.001 ( 0.202 ) 

0.946 ( 

; 0.121) 





F = 1 - 1 

0.188 

0.044 

3.823 ( 0.476 ) 

0.767 ( 

; 0.087 ) 




hn 13 c 


0.479 

0.065 

1.755 ( 0.213 ) 

0.895 ( 

; 0.088 ) 


0.02 

19.74 

h 13 cn 

F = 0-1 

0.199 

0.068 

0.776 ( 0.620 ) 

0.165 ( 

; 0.079 ) 

2.4 

0.07 

26.27 


F = 2-1 

0.479 

0.068 

1.901 ( 0.356 ) 

0.969 ( 

; o.ii8 ) 





F = 1-1 

0.259 

0.068 

1.240 ( 0.409 ) 

0.342 ( 

; 0.083 ) 




hn 13 c 


0.485 

0.059 

1.638 ( 0.221 ) 

0.845 ( 

; 0.083 ) 


0.07 

15.24 

h 13 cn 

F = 0 — 1 

0.130 

0.077 

6.088 ( 1.570 ) 

0.840 ( 

; 0.257 ) 

0.6 

0.10 

14.32 


F = 2-1 

0.317 

0.077 

2.602 ( 0.435 ) 

0.879 ( 

; o.i9o ) 





F = 1-1 

0.177 

0.077 

14.510 ( 3.070 ) 

2.734 ( 

; 0.488 ) 




hn 13 c 


0.629 

0.058 

2.162 ( 0.155 ) 

1.448 ( 

; 0.086 ) 


0.18 

11.34 

h 13 cn 

F = 0-1 

0.461 

0.099 

3.571 ( 0.583 ) 

1.750 ( 

; 0.210) 

32.8 

0.06 

18.62 


F = 2-1 

1.300 

0.099 

3.352 ( 0.250 ) 

4.647 ( 

; 0.429 ) 





F = 1-1 

0.778 

0.099 

5.264 ( 0.578 ) 

4.361 ( 

; 0.453 ) 




hn 13 c 


0.821 

0.072 

3.478 ( 0.168 ) 

3.041 ( 

; o.i3o ) 


0.24 

16.80 

h 13 cn 

F = 0-1 

0.160 

0.068 

1.250 ( 0.384 ) 

0.213 ( 

; o.o7i) 

21.4 

0.16 

18.48 


F = 2-1 

0.405 

0.068 

2.929 ( 0.549 ) 

1.262 ( 

: o.i5i) 





F = 1-1 

0.348 

0.068 

1.341 ( 0.367 ) 

0.496 ( 

; 0.093 ) 




hn 13 c 


0.380 

0.049 

2.127 ( 0.201 ) 

0.861 ( 

; o.o7o ) 


0.06 

16.42 

h 13 cn 

F = 0-1 

0.291 

0.104 

2.932 ( 0.975 ) 

0.908 ( 

; 0.213 ) 

47.6 

0.07 

41.07 


F = 2-1 

0.678 

0.104 

3.728 ( 0.536 ) 

2.691 ( 

; 0.348 ) 





F = 1-1 

0.479 

0.104 

4.074 ( 0.588 ) 

2.078 ( 

; 0.316 ) 




hn 13 c 


0.618 

0.058 

3.305 ( 0.184 ) 

2.173 ( 

; 0.104) 


0.08 

28.39 

h 13 cn 

F = 0-1 

0.207 

0.080 

0.585 ( 0.258 ) 

0.129 ( 

; 0.060 ) 

33.4 

0.06 

47.44 


F = 2 - 1 

0.561 

0.080 

2.265 ( 0.246 ) 

1.352 ( 

; 0.123 ) 





F = 1 - 1 

0.211 

0.080 

3.272 ( 0.704 ) 

0.736 ( 

: o.i4i) 




hn 13 c 


0.430 

0.051 

1.916 ( 0.183 ) 

0.876 ( 

; o.o7o ) 


0.04 

22.45 

h 13 cn 

F = 0 — 1 

0.282 

0.084 

2.651 ( 0.572 ) 

0.796 ( 

: o.i4i) 

59.7 

0.22 

36.16 


F = 2 - 1 

0.891 

0.084 

3.401 ( 0.274 ) 

3.225 ( 

; 0.229 ) 





F = 1 - 1 

0.637 

0.084 

3.706 ( 0.438 ) 

2.512 ( 

; 0.243 ) 




hn 13 c 


0.411 

0.057 

3.157 ( 0.266 ) 

1.382 ( 

; o.ioi) 


0.10 

16.52 

h 13 cn 

F = 0 - 1 

0.146 

0.070 

0.677 ( 0.409 ) 

0.105 ( 

; 0.057 ) 

45.3 

0.07 

31.70 


F = 2 - 1 

0.467 

0.070 

1.471 ( 0.256 ) 

0.731 ( 

; 0.127 ) 





F = 1 - 1 

0.195 

0.070 

9.622 ( 1.590 ) 

1.994 ( 

; 0.270 ) 




hn 13 c 


0.554 

0.060 

1.881 ( 0.153 ) 

1.108 ( 

; 0.080 ) 


0.05 

19.83 

h 13 cn 

F = 0 - 1 

0.141 

0.067 

7.066 ( 1.720 ) 

1.059 ( 

; 0.235 ) 

43.8 

0.11 

71.79 


F = 2 - 1 

1.140 

0.067 

2.667 ( 0.173 ) 

3.233 ( 

; 0.213 ) 





F = 1 - 1 

0.669 

0.067 

3.789 ( 0.266 ) 

2.699 ( 

; 0.155 ) 




hn 13 c 


0.672 

0.072 

2.401 ( 0.179 ) 

1.718 ( 

; o.no ) 


0.07 

29.08 

h 13 cn 

F = 0 - 1 

0.192 

0.067 

0.434 ( 4.430 ) 

0.088 ( 

; 0.042 ) 

76.5 

0.10 

23.04 


F = 2 - 1 

0.440 

0.067 

2.458 ( 0.296 ) 

1.152 ( 

; 0.133 ) 





IRAS 18308-0841 





46 


Tabic 8—Continued 


Source Name 

Molecule 

Transition 

Tmb 

Trms 

Av 


^LSR 

Tthin 

Tex 




[K] 

[K] 

[kms -1 ] 

[K kms -1 ] 

[kms -1 ] 


[K] 




F = 1-1 

0.194 

0.067 

3.688 ( 0.800 ) 

0.760 ( 0.149 ) 





hn 13 c 


0.508 

0.049 

2.366 ( 0.165 ) 

1.280 ( 0.075 ) 


0.14 

11.67 

IRAS 18345-0641 

h 13 cn 

F = 0-1 

0.104 

0.060 

6.129 ( 2.100 ) 

0.678 ( 0.216 ) 

95.5 

0.06 

21.09 



F = 2-1 

0.329 

0.060 

3.105 ( 0.568 ) 

1.087 ( 0.245 ) 






F = 1-1 

0.238 

0.060 

4.701 ( 1.170 ) 

1.189 ( 0.217 ) 





hn 13 c 


0.522 

0.050 

1.859 ( 0.149 ) 

1.032 ( 0.069 ) 


0.07 

17.37 

IRAS 18440-0148 

h 13 cn 

F = 0-1 

0.093 

0.067 

0.955 ( 0.533 ) 

0.095 ( 0.058 ) 

97.7 

0.08 

16.09 



F = 2-1 

0.292 

0.067 

2.385 ( 0.349 ) 

0.741 ( 0.100 ) 






F = 1-1 

0.120 

0.067 

2.109 ( 1.010 ) 

0.270 ( 0.101 ) 





hn 13 c 


0.278 

0.062 

1.414 ( 0.382 ) 

0.419 ( 0.082 ) 


0.03 

6.15 

IRAS 18470-0044 

h 13 cn 

F = 0-1 

0.062 

0.075 

6.469 ( 3.490 ) 

0.424 ( 0.195 ) 

96.3 

0.16 

21.42 



F = 2 - 1 

0.467 

0.075 

3.026 ( 0.378 ) 

1.505 ( 0.158 ) 






F = 1 - 1 

0.253 

0.075 

2.999 ( 0.605 ) 

0.807 ( 0.137 ) 





hn 13 c 


0.266 

0.060 

2.588 ( 0.459 ) 

0.733 ( 0.101 ) 


0.06 

11.77 

IRAS 18511+0146 

h 13 cn 

F = 0 - 1 

0.156 

0.061 

0.586 ( 0.365 ) 

0.098 ( 0.047 ) 

57.0 

0.09 

9.46 



F = 2 - 1 

0.373 

0.061 

2.540 ( 0.405 ) 

1.010 ( 0.232 ) 






F = 1 - 1 

0.168 

0.061 

7.099 ( 1.830 ) 

1.267 ( 0.299 ) 





hn 13 c 


0.759 

0.076 

1.826 ( 0.150 ) 

1.476 ( 0.103 ) 


0.12 

15.23 

IRAS 18530+0215 

h 13 cn 

F = 0 - 1 

0.166 

0.076 

2.664 ( 0.636 ) 

0.472 ( 0.115 ) 

77.3 

0.10 

33.43 



F = 2 - 1 

0.621 

0.076 

2.694 ( 0.251 ) 

1.782 ( 0.134 ) 






F = 1 - 1 

0.431 

0.076 

2.908 ( 0.356 ) 

1.333 ( 0.136 ) 





hn 13 c 


0.560 

0.071 

2.727 ( 0.207 ) 

1.626 ( 0.113 ) 


0.09 

21.19 

IRAS 19220+1432 

h 13 cn 

F = 0 - 1 





69.8 

0.02 

23.75 


F = 2-1 
F = 1 - 1 


IRAS 19410+2336 


IRAS 19411+2306 


IRAS20126+4104 


hn 13 c 



0.234 a 

h 13 cn 

F = 0 - 

1 

0.138 


F = 2 - 

1 

1.240 


F = 1 - 

1 

0.547 

hn 13 c 



0.903 

h 13 cn 

F = 0 - 

1 

0.086 


F = 2 - 

1 

0.366 


F = 1 - 

1 

0.138 

hn 13 c 



0.335 

h 13 cn 

F = 0 - 

1 

0.138 


F = 2 - 

1 

0.639 


F = 1 - 

1 

0.387 

hn 13 c 



1.380 

h 13 cn 

F = 0 - 

1 

0.132 


F = 2- 

1 

0.315 


F = 1 - 

1 

0.212 

hn 13 c 



0.462 


0.081 

3.164 ( 0.503 ) 

0.787 ( 

0.154 

3.219 ( 3.590 ) 

0.474 ( 

0.154 

2.660 ( 0.275 ) 

3.502 ( 

0.154 

3.276 ( 0.682 ) 

1.908 ( 

0.063 

1.785 ( 0.097 ) 

1.715 ( 

0.061 

1.457 ( 0.467 ) 

0.134 ( 

0.061 

1.717 ( 0.233 ) 

0.669 ( 

0.061 

4.085 ( 1.010 ) 

0.600 ( 

0.066 

1.127 ( 0.249 ) 

0.401 ( 

0.056 

6.360 ( 1.270 ) 

0.931 ( 

0.056 

2.626 ( 0.204 ) 

1.786 ( 

0.056 

2.926 ( 0.281 ) 

1.206 ( 

0.098 

1.989 ( 0.106 ) 

2.928 ( 

0.047 

0.617 ( 0.387 ) 

0.087 ( 

0.047 

2.396 ( 0.238 ) 

0.804 ( 

0.047 

1.921 ( 0.385 ) 

0.433 ( 

0.078 

2.261 ( 0.269 ) 

1.113 ( 


0.132 ) 


0.08 

9.26 

0.344 ) 

22.6 

0.12 

53.47 

0.290 ) 
0.316 ) 
0.082 ) 


0.09 

23.21 

0.064 ) 

29.2 

0.05 

28.96 

0.081 ) 
0.124 ) 
0.072 ) 


0.05 

12.39 

0.157 ) 

-3.9 

0.07 

58.92 

0.121 ) 
0.097 ) 
0.136 ) 


0.12 

26.63 

0.038 ) 

11.7 

0.03 

43.72 

0.072 ) 
0.068 ) 
0.116 ) 


0.04 

15.71 


IRAS20343+4129 
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Table 8—Continued 


Source Name 

Molecule 

Transition 

T MB 

[K] 

Irms 

[K] 

Av 

[kms -1 ] 

/ TmB dv 
[K kms -1 ] 

V LSR 

[kms -1 ] 

T thin 

Tex 

[K] 

IRAS22198+6336 

h 13 cn 

F 

= 0 - 

1 

0.214 

0.048 

0.660 ( 0.203 ) 

0.150 ( 0.040 ) 

-11.0 

0.10 

14.12 



F 

= 2- 

1 

0.309 

0.048 

1.400 ( 0.253 ) 

0.461 ( 0.077 ) 






F 

= 1 - 

1 

0.095 

0.048 

4.848 ( 1.480 ) 

0.488 ( 0.118 ) 





hn 13 c 




0.326 

0.041 

1.302 ( 0.165 ) 

0.451 ( 0.047 ) 


0.05 

9.92 

IRAS23033+5951 

h 13 cn 

F 

= 0 - 

1 

0.206 

0.063 

1.716 ( 0.371 ) 

0.377 ( 0.078 ) 

-53.7 

0.10 

33.53 



F 

= 2- 

1 

0.566 

0.063 

2.796 ( 0.219 ) 

1.685 ( 0.113 ) 






F 

= 1 - 

1 

0.305 

0.063 

3.277 ( 0.471 ) 

1.063 ( 0.125 ) 





hn 13 c 




0.283 

0.060 

3.246 ( 0.483 ) 

0.976 ( 0.113 ) 


0.06 

14.46 


— “The detection of the HN 13 C 


line can be confirmed by eye even though their SNR is smaller than 3. 


Note. 
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Table 9. UCHII Observed Line Parameters 


Source Name 

Molecule 

Transition 

%B 

Irms 

Av 

/ TMBdv 

^LSR 

T thin 

Tex 




[K] 

[K] 

[kms" 1 ] 

[K kms -1 ] 

[kms -1 ] 


[K] 


IRAS02232+6138 


IRAS02575+6017 


IRAS03035+5819 


IRAS06053-0622 


IRAS06084-0611 


IRAS06056+2131 


IRAS06058+2138 


IRAS06061+2151 


IRAS06099+1800 


IRAS 17574-2403 


h 13 cn 

F = 0 - 1 

0.288 

0.042 

3.422 1 

[ 0.361 ) 

1.049 ( 0.083 ) 

-46.4 

0.15 

59.33 


F = 2-1 

1.160 

0.042 

3.454 1 

[ 0.093 ) 

4.269 ( 0.096 ) 





F = 1 - 1 

0.658 

0.042 

3.381 I 

{ 0.164 ) 

2.369 ( 0.095 ) 




hn 13 c 


0.378 

0.039 

2.938 I 

{ 0.185 ) 

1.182 ( 0.065 ) 


0.05 

29.14 

h 13 cn 

F = 0 - 1 

0.096 

0.048 

2.062 i 

{ 0.662 ) 

0.211 ( 0.065 ) 

-38.0 

0.04 

47.81 


F = 2-1 

0.543 

0.048 

2.347 I 

{ 0.154 ) 

1.356 ( 0.074 ) 





F = 1 - 1 

0.228 

0.048 

2.396 I 

{ 0.392 ) 

0.583 ( 0.076 ) 




hn 13 c 


0.300 

0.041 

1.381 i 

{ 0.154 ) 

0.441 ( 0.046 ) 


0.01 

15.98 

h 13 cn 

F = 0-1 

0.111 

0.044 

1.035 I 

{ 0.543 ) 

0.122 ( 0.048 ) 

-39.3 

0.06 

25.20 


F = 2-1 

0.403 

0.044 

1.619 i 

{ 0.156 ) 

0.695 ( 0.057 ) 





F = 1 - 1 

0.250 

0.044 

1.826 i 

{ 0.316 ) 

0.487 ( 0.064 ) 




hn 13 c 


0.167 

0.040 

1.086 i 

{ 0.263 ) 

0.193 ( 0.041 ) 


0.04 

9.39 

h 13 cn 

F = 0 - 1 

0.165 

0.046 

0.900 I 

{ 0.336 ) 

0.158 ( 0.047 ) 

10.4 

0.06 

35.59 


F = 2-1 

0.352 

0.046 

3.324 i 

{ 0.253 ) 

1.245 ( 0.084 ) 





F = 1 - 1 

0.269 

0.046 

1.849 i 

{ 0.230 ) 

0.530 ( 0.062 ) 




hn 13 c 


0.141“ 

0.051 

1.530 I 

' 0.373 ) 

0.229 ( 0.057 ) 


0.02 

16.41 

h 13 cn 

F = 0-1 

0.181 

0.046 

5.264 i 

; 1.070 ) 

1.015 ( 0.136 ) 

11.3 

0.11 

66.14 


F = 2-1 

1.110 

0.046 

2.914 i 

; 0.095 ) 

3.455 ( 0.091 ) 





F = 1 - 1 

0.646 

0.046 

2.719 l 

{ 0.117 ) 

1.871 ( 0.076 ) 




hn 13 c 


0.374 

0.043 

2.294 1 

{ 0.191 ) 

0.913 ( 0.065 ) 


0.04 

24.72 

h 13 cn 

F = 0 — 1 

0.177 

0.049 

2.173 1 

{ 0.374 ) 

0.408 ( 0.068 ) 

2.8 

0.06 

49.07 


F = 2-1 

0.615 

0.049 

2.281 1 

( 0.137 ) 

1.493 ( 0.076 ) 





F = 1-1 

0.329 

0.049 

2.543 l 

{ 0.381 ) 

0.892 ( 0.093 ) 




hn 13 c 


0.359 

0.051 

2.594 1 

' 0.273 ) 

0.991 ( 0.085 ) 


0.06 

21.60 

h 13 cn 

F = 0 - 1 

0.197 

0.046 

2.558 l 

{ 0.398 ) 

0.535 ( 0.073 ) 

3.4 

0.07 

54.80 


F = 2-1 

0.641 

0.046 

2.733 I 

{ 0.139 ) 

1.864 ( 0.081 ) 





F = 1-1 

0.341 

0.046 

2.995 I 

{ 0.317 ) 

1.087 ( 0.090 ) 




hn 13 c 


0.409 

0.047 

2.054 i 

{ 0.181 ) 

0.895 ( 0.067 ) 


0.05 

22.87 

h 13 cn 

F = 0 — 1 

0.134 

0.038 

1.605 I 

{ 0.250 ) 

0.230 ( 0.043 ) 

-0.8 

0.09 

23.62 


F = 2-1 

0.374 

0.038 

2.568 I 

{ 0.182 ) 

1.021 ( 0.061 ) 





F = 1-1 

0.267 

0.038 

2.795 I 

{ 0.249 ) 

0.796 ( 0.062 ) 




hn 13 c 


0.130 

0.042 

1.590 I 

{ 0.352 ) 

0.220 ( 0.050 ) 


0.03 

10.51 

h 13 cn 

F = 0 - 1 

0.364 

0.056 

1.214 i 

{ 0.207 ) 

0.470 ( 0.065 ) 

7.3 

0.10 

72.67 


F = 2 - 1 

1.410 

0.056 

2.348 l 

{ 0.084 ) 

3.534 ( 0.100 ) 





F = 1 - 1 

0.655 

0.056 

3.473 I 

{ 0.248 ) 

2.422 ( 0.126 ) 




hn 13 c 


0.340 

0.045 

1.834 i 

{ 0.213 ) 

0.664 ( 0.062 ) 


0.04 

20.68 

h 13 cn 

F = 0 - 1 

0.807 

0.074 

3.990 l 

' 0.301 ) 

3.427 ( 0.270 ) 

8.9 

0.66 

77.78 


F = 2 - 1 

2.240 

0.074 

2.588 I 

' 0.105 ) 

6.157 ( 0.435 ) 





F = 1 - 1 

2.570 

0.074 

9.511 i 

{ 0.226 ) 

26.070 ( 0.723 ) 




hn 13 c 


2.120 

0.078 

3.813 I 

{ 0.090 ) 

8.624 ( 0.158 ) 


0.12 

87.26 

h 13 cn 

F = 0 - 1 

0.116 

0.072 

4.091 I 

' 1.540 ) 

0.503 ( 0.161 ) 

18.9 

0.15 

25.85 


IRAS 17599-2148 





49 


Tabic 9—Continued 


Source Name 

Molecule 

Transition 

t mb 

[K] 

^rms 

[K] 

Av 

[kms -1 ] 

f %B dv 
[K kms -1 ] 

^LSR 

[kms -1 ] 

T thin 

^ex 

[K] 



F = 2-1 

0.359 

0.072 

6.443 ( 0.845 ) 

2.461 ( 0.255 ) 






F = 1-1 

0.225 

0.072 

2.042 ( 0.641 ) 

0.489 ( 0.162 ) 





hn 13 c 


0.465 

0.053 

4.160 ( 0.242 ) 

2.057 ( 0.105 ) 


0.14 

19.96 

IRAS 18032-2137 

h 13 cn 

F = 0-1 

0.402 

0.073 

4.532 ( 0.480 ) 

1.940 ( 0.169 ) 

35.1 

0.48 

35.02 



F = 2-1 

1.520 

0.073 

3.394 ( 0.148 ) 

5.479 ( 0.238 ) 






F = 1-1 

0.960 

0.073 

4.667 ( 0.349 ) 

4.769 ( 0.270 ) 





hn 13 c 


1.520 

0.065 

4.450 ( 0.111 ) 

7.183 ( 0.141 ) 


0.24 

37.29 

IRAS 18032-2032 

h 13 cn 

F = 0 - 1 







3.8 

0.09 

43.34 



F = 2 - 1 












F = 1 - 1 











hn 13 c 


0.179 

0.053 

4.476 ( 0.690 ) 

0.852 ( 0.111 ) 


0.03 

33.77 

IRAS 18075-1956 

h 13 cn 

F = 0 - 1 







-2.0 

0.04 

35.36 



F = 2 - 1 












F = 1 - 1 











hn 13 c 


0.206 

0.053 

2.869 

( 0.544 ) 

0.628 ( 

; 0.094 ) 


0.04 

23.85 

IRAS 18162-2048 

h 13 cn 

F = 0 - 1 

0.193 

0.050 

1.607 

( 0.325 ) 

0.331 ( 

; 0.061 ) 

12.5 

0.15 

59.33 



F = 2 - 1 

0.666 

0.050 

2.799 

( 0.140 ) 

1.984 ( 

; 0.084 ) 






F = 1 - 1 

0.425 

0.050 

2.281 

( 0.207 ) 

1.033 ( 

; 0.077 ) 





hn 13 c 


0.241 

0.052 

2.485 

( 0.316 ) 

0.638 ( 

; 0.077 ) 


0.05 

14.75 

IRAS 18100-1854 

h 13 cn 

F = 0 - 1 

0.242 

0.048 

5.892 

( 0.927 ) 

1.515 ( 

; 0.232 ) 

38.4 

0.04 

47.81 



F = 2 - 1 

0.546 

0.048 

3.824 

( 0.276 ) 

2.222 ( 

; 0.199 ) 






F = 1 - 1 

1.050 

0.048 

4.710 

( 0.307 ) 

5.251 ( 

; 0.339 ) 





hn 13 c 


1.010 

0.056 

3.972 

( 0.109 ) 

4.281 ( 

; 0.106 ) 


0.21 

23.27 

IRAS 18174-1612 

h 13 cn 

F = 0 - 1 

0.213 

0.051 

2.874 

( 0.418 ) 

0.651 ( 

; 0.086 ) 

19.4 

0.06 

25.20 



F = 2 - 1 

0.869 

0.051 

3.121 

( 0.144 ) 

2.887 ( 

; o.io7) 






F = 1 - 1 

0.502 

0.051 

3.096 

( 0.232 ) 

1.656 ( 

; o.io4) 





hn 13 c 


0.440 

0.053 

2.949 

( 0.263 ) 

1.380 ( 

; 0.096 ) 




IRAS 18317-0757 

h 13 cn 

F = 0 - 1 

0.156 

0.052 

3.424 

( 0.596 ) 

0.568 ( 

; 0.093 ) 

80.1 

0.23 

19.35 



F = 2 - 1 

0.697 

0.052 

2.936 

( 0.162 ) 

2.177 ( 

; 0.099 ) 






F = 1 - 1 

0.291 

0.052 

3.473 

( 0.376 ) 

1.077 ( 

; o.ioi ) 





hn 13 c 


0.245 

0.043 

2.859 

( 0.285 ) 

0.744 ( 

; 0.069 ) 




IRAS 18316-0602 

h 13 cn 

F = 0 - 1 

0.227 

0.048 

5.739 

( 1.050 ) 

1.389 ( 

; 0.223 ) 

42.4 

0.14 

65.58 



F = 2 - 1 

0.509 

0.048 

4.436 

( 0.370 ) 

2.402 ( 

; 0.495 ) 






F = 1 - 1 

0.536 

0.048 

7.546 

( 0.774 ) 

4.308 ( 

; 0.550 ) 





hn 13 c 


0.456 

0.049 

3.717 

( 0.227 ) 

1.805 ( 

; 0.094 ) 




IRAS 18403-0417 

h 13 cn 

F = 0 - 1 







95.7 

0.33 

43.01 



F = 2 - 1 












F = 1 - 1 











hn 13 c 


0.711 

0.051 

4.521 ( 0.183 ) 

3.423 ( 0.111 ) 


0.13 

32.96 


h 13 cn F = 0 - 1 
F = 2-1 
F = 1-1 


IRAS 18479-0005 


17.5 


0.16 


40.30 
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Table 9—Continued 


Source Name 

Molecule 

Transition 

^MB 

[K] 

Trms 

[K] 

Av 

[kms -1 ] 

/ iMBdv 
[K kms -1 ] 

^LSR 

[kms -1 ] 

T thin 

Tex 

[K] 


hn 13 c 


0.196 

0.065 

4.882 ( 1.640 ) 

1.019 ( 0.207 ) 


0.04 

26.25 

IRAS 18507+0110 

h 13 cn 

F = 0 - 1 

0.913 

0.086 

6.074 ( 0.344 ) 

5.901 ( 0.263 ) 

58.0 

1.58 

40.96 



F = 2-1 

3.400 

0.086 

4.625 ( 0.061 ) 

16.750 ( 0.211 ) 






F = 1 - 1 

1.990 

0.086 

3.809 ( 0.126 ) 

8.071 ( 0.239 ) 





hn 13 c 


0.930 

0.060 

5.500 ( 0.161 ) 

5.444 ( 0.138 ) 


0.11 

51.10 

IRAS 19095+0930 

h 13 cn 

F = 0-1 






44.1 

0.14 

37.08 



F = 2-1 











F = 1 — 1 










hn 13 c 


0.212 

0.056 

4.875 

( 0.688 ) 

1.100 ( 0.126 ) 


0.06 

21.97 

IRAS 19110+1045 

h 13 cn 

F = 0 — 1 

0.185 

0.053 

4.615 

( 0.434 ) 

0.910 ( 0.081 ) 

59.5 

0.15 

41.90 



F = 2 — 1 

0.528 

0.053 

4.406 

( 0.434 ) 

2.478 ( 0.081 ) 






F = 1 — 1 

0.447 

0.053 

4.744 

( 0.434 ) 

2.258 ( 0.081 ) 





hn 13 c 


0.234 

0.058 

2.246 

( 0.501 ) 

0.559 ( 0.094 ) 


0.02 

20.77 

IRAS20081+3122 

h 13 cn 

F = 0 — 1 

0.166 

0.041 

3.794 

( 0.495 ) 

0.669 ( 0.079 ) 

11.7 

0.16 

34.28 



F = 2-1 

0.642 

0.041 

3.823 

( 0.210 ) 

2.611 ( 0.126 ) 






F = 1-1 

0.372 

0.041 

4.018 

( 0.388 ) 

1.590 ( 0.129 ) 





hn 13 c 


0.738 

0.060 

3.936 

( 0.182 ) 

3.094 ( 0.120 ) 


0.13 

26.94 

IRAS20255+3712 

h 13 cn 

F = 0 — 1 

0.096 

0.051 

1.692 

( 0.646 ) 

0.173 ( 0.069 ) 

-1.5 

0.04 

47.82 



F = 2-1 

0.403 

0.051 

2.142 

( 0.253 ) 

0.920 ( 0.082 ) 






F = 1 — 1 

0.243 

0.051 

2.734 

( 0.503 ) 

0.708 ( 0.095 ) 





hn 13 c 


0.440 

0.048 

2.003 

( 0.211 ) 

0.939 ( 0.074 ) 


0.09 

14.77 

IRAS20178+4046 

h 13 cn 

F = 0 - 1 

0.072 

0.044 

2.367 

( 0.673 ) 

0.183 ( 0.061 ) 

1.1 

0.07 

22.42 



F = 2 - 1 

0.302 

0.044 

1.601 

( 0.186 ) 

0.515 ( 0.054 ) 






F = 1 - 1 

0.130 

0.044 

2.489 

( 0.447 ) 

0.345 ( 0.065 ) 





hn 13 c 


0.347 

0.060 

1.886 

( 0.247 ) 

0.697 ( 0.081 ) 


0.07 

10.44 

IRAS20350+4126 

h 13 cn 

F = 0 - 1 

0.199 

0.052 

1.111 

( 0.862 ) 

0.236 ( 0.104 ) 

-2.5 

0.15 

20.18 



F = 2 - 1 

0.480 

0.052 

2.453 

( 0.219 ) 

1.254 ( 0.087 ) 






F = 1 - 1 

0.320 

0.052 

2.509 

( 0.326 ) 

0.853 ( 0.088 ) 





hn 13 c 


0.349 

0.051 

2.357 

( 0.267 ) 

0.876 ( 0.081 ) 


0.09 

12.14 

IRAS22543+6145 

h 13 cn 

F = 0 — 1 

0.208 

0.052 

2.092 

( 0.379 ) 

0.463 ( 0.075 ) 

-10.6 

0.08 

70.60 



F = 2 - 1 

0.780 

0.052 

3.101 

( 0.144 ) 

2.573 ( 0.099 ) 






F = 1 - 1 

0.469 

0.052 

3.024 

( 0.246 ) 

1.511 ( 0.099 ) 





hn 13 c 


0.548 

0.082 

2.762 

( 0.277 ) 

1.611 ( 0.136 ) 


0.04 

27.72 

IRAS23138+5945 

h 13 cn 

F = 0 - 1 

0.119 

0.050 

1.039 

( 0.478 ) 

0.132 ( 0.052 ) 

-44.4 

0.04 

45.71 



F = 2 - 1 

0.295 

0.050 

3.280 

( 0.317 ) 

1.031 ( 0.089 ) 






F = 1 - 1 

0.281 

0.050 

2.072 

( 0.288 ) 

0.621 ( 0.073 ) 





hn 13 c 


0.250 

0.062 

2.094 

( 0.395 ) 

0.557 ( 0.089 ) 


0.04 

16.76 

IRAS23133+6050 

h 13 cn 

F = 0 - 1 

0.157 

0.052 

1.907 

( 0.545 ) 

0.319 ( 0.073 ) 

-56.7 

0.07 

48.41 



F = 2 - 1 

0.718 

0.052 

2.557 

( 0.140 ) 

1.954 ( 0.086 ) 






F = 1 - 1 

0.377 

0.052 

1.970 

( 0.253 ) 

0.790 ( 0.078 ) 





hn 13 c 


0.150“ 

0.062 

3.928 

( 0.883 ) 

0.626 ( 0.121 ) 


0.01 

17.34 






Note. 


The detection of the HN 13 C line can be confirmed by eye even though their SNR is smaller than 3. 
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Table 10. Derived mean values and standard deviations of r t hi n and T ex for each 

evolutionary phases. 


Evolutionary Phases 

t h 13 cn 

Tex(HCN) 

[K] 

t hn 13 c 

Tex(HNC) 

[K] 

qIRDCc 

0.11 (0.03) 

18.79 (8.74) 

0.11 (0.03) 

15.10 (2.58) 

aIRDCc 

0.16 (0.08) 

22.38 (13.33) 

0.13 (0.08) 

15.77 (7.03) 

HMPO 

0.09 (0.05) 

30.88 (15.77) 

0.08 (0.05) 

16.99 (6.11) 

UCHII 

0.18 (0.29) 

43.92 (16.12) 

0.07 (0.06) 

24.67 (15.26) 


Table 11. 

QIRDC Derived Line Parameters 


Source Name 

N(HCN) 

[10 14 cm -2 ] 

N(HNC) 

[10 14 cm -2 ] 

HCN 

HNC 

N(H 2 ) 

[10 22 cm“ 2 ] 

X(HCN) 

[io- 8 ] 

X(HNC) 

[IO" 8 ] 

G018.82.MM6 

i.i ( 

0.11 ) 

i.i ( 

O 

o 

CO 

1.06 ( 0.156 ) 




G028.37.MM9 

1.2 ( 

i> 

o 

o 

1.4 ( 

0.05 ) 

0.87 ( 0.069 ) 

1.75 

0.7 

0.8 

G028.37.MM17 

1.5 ( 

0.12 ) 

1.5 ( 

t- 

o 

o 

1.02 ( 0.141 ) 

0.68 

2.2 

2.1 

G028.53.MM7 

0.8 ( 

0.13 ) 

1.0 ( 

0.10 ) 

0.81 ( 0.130 ) 




G031.97.MM9 

2.0 ( 

0.07 ) 

1.5 ( 

0.05 ) 

1.30 ( 0.116 ) 

1.27 

1.5 

1.2 

G033.69.MM8 

1.0 ( 

0.12 ) 

0.9 ( 

0.07 ) 

1.17 ( 0.162 ) 

0.55 

1.9 

1.6 

G034.77.MM4 

1.0 ( 

0.16 ) 

1.3 ( 

0.06 ) 

0.77 ( 0.132 ) 




G035.39.MM5 

1.1 ( 

0.10 ) 

1.5 ( 

0.06 ) 

0.75 ( 0.087 ) 





Table 12. 

AIRDC Derived Line Parameters 


Source Name 

N(HCN) 

[10 14 cm -2 ] 

N(HNC) 

[10 14 cm -2 ] 


HCN 

HNC 

N(H 2 ) 

[10 22 cm -2 ] 

X(HCN) 

[IO" 8 ] 

X(HNC) 

[IO" 8 ] 

G019.27.MM2 

1.3 ( 

; 0.08 ) 

0.7 ( 

0.08 ) 

1.90 

( 0.210 ) 

0.59 

2.2 

1.1 

G022.35.MM1 

2.7 ( 

; 0.06 ) 

1.0 ( 

0.07 ) 

2.75 

( 0.245 ) 




G023.60.MM1 

6.8 ( 

; o.o3 ) 

2.1 ( 

0.03 ) 

3.24 

( 0.127 ) 




G024.60.MM1 

2.8 ( 

; o.o6 ) 

0.7 ( 

0.10 ) 

4.20 

( 0.492 ) 

1.08 

2.6 

0.6 

G028.37.MM4 

3.8 ( 

: o.o4 ) 

2.9 ( 

0.04 ) 

1.31 

( 0.065 ) 

1.53 

2.5 

1.9 

G028.37.MM6 

2.5 ( 

: o.o4 > 

1.2 ( 

0.04 ) 

2.01 

( 0.115 ) 

1.13 

2.2 

1.1 

G030.97.MM1 

3.2 ( 

: o.o5 > 

1.3 ( 

0.05 ) 

2.34 

( 0.173 ) 

1.04 

3.0 

1.3 

G031.97.MM1 

23.1 ( 

: o.oi > 

3.6 ( 

0.03 ) 

6.45 

( 0.187 ) 

3.51 

6.6 

1.0 

G033.69.MM4 

1.8 ( 

: o.o7 > 

1.6 ( 

0.03 ) 

1.12 

( 0.085 ) 

1.08 

1.6 

1.5 

G033.69.MM5 

0.6 ( 

; 0.09 ) 

0.8 ( 

0.07 ) 

0.81 

( 0.094 ) 




G034.43.MM1 

27.2 ( 

: o.o2 > 

4.2 ( 

0.03 ) 

6.47 

( 0.230 ) 

5.55 

4.9 

0.8 

G034.43.MM3 

8.9 ( 

; 0.03 ) 

3.0 ( 

0.03 ) 

2.93 

( 0.117 ) 

1.18 

7.6 

2.6 

G034.43.MM4 

10.0 ( 

: o.o2) 

3.5 ( 

0.03 ) 

2.85 

( 0.084 ) 

3.16 

3.2 

1.1 

G035.39.MM7 

0.9 ( 

: o.os) 

1.4 ( 

0.05 ) 

0.63 

( 0.061 ) 




G048.65.MM1 

0.6 ( 

: o.i5) 

0.4 ( 

0.15 ) 

1.65 

( 0.344 ) 

0.33 

1.8 

1.1 

G053.25.MM6 

0.7 ( 

: o.i5) 

0.4 ( 

0.18 ) 

1.76 

( 0.401 ) 
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Table 13. HMPO Derived Line Parameters 


Source Name 

N(HCN) 

[10 14 cm -2 ] 

N(HNC) 

[10 14 cm -2 ] 


HCN 

HNC 

N(H 2 ) 

[10 22 cm~ 2 ] 

X(HCN) 

[lO" 8 ] 

X(HNC) 

[lO" 8 ] 

IRAS05358+3543 

11.2 ( 

: o.o4) 

3.5 ( 

0.11 ) 

3.22 

( 0.376 ) 

2.54 

4.4 

1.4 

IRAS05373+2349 

4.5 ( 

: o.o7) 

1.8 ( 

0.09 ) 

2.44 

( 0.276 ) 

0.85 

5.3 

2.2 

IRAS 18024-2119 

5.6 ( 

: o.o4) 

1.9 ( 

0.05 ) 

2.94 

( 0.186 ) 

1.44 

3.9 

1.3 

IRAS 18089-1732 

13.4 ( 

: o.o2) 

3.2 ( 

0.04 ) 

4.25 

( 0.193 ) 

3.29 

4.1 

1.0 

IRAS 18102-1800 

3.1 ( 

; 0.05 ) 

1.1 ( 

0.06 ) 

2.95 

( 0.234 ) 

2.00 

1.6 

0.5 

IRAS 18144-1723 

10.9 ( 

: o.o4) 

3.2 ( 

0.04 ) 

3.44 

( 0.186 ) 

1.13 

9.7 

2.8 

IRAS 18151-1208 

6.9 ( 

; 0.08 ) 

1.3 ( 

0.07 ) 

5.21 

( 0.536 ) 

1.67 

4.2 

0.8 

IRAS 18182-1433 

12.2 ( 

; 0.03 ) 

1.4 ( 

0.06 ) 

8.98 

( 0.604 ) 

2.00 

6.1 

0.7 

IRAS 18223-1243 

4.8 ( 

; 0.05 ) 

1.2 ( 

0.08 ) 

3.95 

( 0.367 ) 

1.41 

3.4 

0.9 

IRAS 18264-1152 

27.8 ( 

: o.o2) 

3.0 ( 

0.06 ) 

9.27 

( 0.612 ) 

2.29 

12.2 

1.3 

IRAS 18308-0841 

2.8 ( 

: o.o6) 

0.9 ( 

0.06 ) 

3.02 

( 0.255 ) 

0.92 

3.1 

1.0 

IRAS 18345-0641 

3.0 ( 

: o.o4) 

1.0 ( 

0.05 ) 

2.83 

( 0.188 ) 

0.68 

4.4 

1.5 

IRAS 18440-0148 

1.0 ( 

: o.n) 

0.3 ( 

0.16 ) 

3.95 

( 0.768 ) 

0.42 

2.4 

0.6 

IRAS 18470-0044 

2.8 ( 

: o.o7) 

0.6 ( 

0.09 ) 

4.67 

( 0.546 ) 

0.98 

2.8 

0.6 

IRAS 18511+0146 

2.3 ( 

: o.o6) 

1.6 ( 

0.06 ) 

1.37 

( 0.110 ) 

1.46 

1.5 

1.1 

IRAS 18530+0215 

7.2 ( 

; 0.05 ) 

2.1 ( 

0.06 ) 

3.41 

( 0.256 ) 

2.68 

2.7 

0.8 

IRAS 19220+1432 

1.5 ( 

: o.i6) 

0.6 ( 

0.15 ) 

2.57 

( 0.571 ) 

1.01 

1.5 

0.6 

IRAS 19410+2336 

23.9 ( 

: o.o5) 

3.1 ( 

0.04 ) 

7.67 

( 0.484 ) 

2.12 

11.3 

1.5 

IRAS 19411+2306 

2.9 ( 

: o.io) 

0.5 ( 

0.15 ) 

5.58 

( 1.009 ) 

0.56 

5.2 

0.9 

IRAS20126+4104 

18.6 ( 

: o.o3) 

6.4 ( 

0.06 ) 

2.89 

( 0.187 ) 

2.07 

9.0 

3.1 

IRAS20343+4129 

4.5 ( 

: o.o7) 

1.5 ( 

0.10 ) 

3.02 

( 0.363 ) 

1.31 

3.4 

1.1 

IRAS22198+6336 

1.7 ( 

; o.o8 ) 

0.6 ( 

0.09 ) 

2.67 

( 0.321 ) 




IRAS23033+5951 

10.8 ( 

: o.o4) 

1.9 ( 

0.07 ) 

5.61 

( 0.463 ) 

1.48 

7.3 

1.3 
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Tabic 14. UCHII Derived Line Parameters 


Source Name 

N(HCN) 

[10 14 cm -2 ] 

N(HNC) 

[10 14 cm -2 ] 

HCN 

HNC 

N(H 2 ) 

[10 22 cm -2 ] 

X(HCN) 

[io- 8 ] 

X(HNC) 

[IO" 8 ] 

IRAS02232+6138 

48.3 

( 0.01 ) 

3.7 ( 

; 0.05 ) 

12.99 ( 

; 0.650 ) 

3.88 

12.4 

1.0 

IRAS02575+6017 

10.8 

( 0.04 ) 

0.8 ( 

: o.i3 ) 

13.83 ( 

; 1.803 ) 

0.57 

19.0 

1.4 

IRAS03035+5819 

4.5 

( 0.04 ) 

0.2 ( 

: 0.19 ) 

18.08 ( 

; 3.542 ) 

0.50 

8.8 

0.5 

IRAS06053-0622 

6.3 

( 0.05 ) 

0.4 ( 

: o.i6 ) 

15.35 ( 

; 2.573 ) 




IRAS06084-0611 

33.0 

( 0.01 ) 

2.8 ( 

; 0.03 ) 

11.88 ( 

; 0.435 ) 

1.80 

18.3 

1.5 

IRAS06056+2131 

11.5 

( 0.04 ) 

2.0 ( 

; 0.06 ) 

5.86 ( 

; 0.386 ) 

1.53 

7.6 

1.3 

IRAS06058+2138 

16.8 

( 0.03 ) 

2.0 ( 

; 0.06 ) 

8.35 ( 

; 0.561 ) 

1.49 

11.3 

1.4 

IRAS06061+2151 

4.6 

( 0.04 ) 

0.2 ( 

: o.i7) 

21.27 ( 

; 3.686 ) 

0.48 

9.7 

0.5 

IRAS06099+1800 

46.1 

( 0.02 ) 

1.6 ( 

: o.o9) 

29.22 ( 

; 2.620 ) 

1.29 

35.8 

1.2 

IRAS 17574-2403 

166.0 

( 0.01 ) 

45.2 ( 

: o.oi) 

3.68 ( 

; 0.056 ) 

3.56 

46.7 

12.7 

IRAS 17599-2148 

5.2 

( 0.04 ) 

2.4 ( 

: o.o4) 

2.15 ( 

; o.ii9) 

1.51 

3.4 

1.6 

IRAS 18032-2137 

21.8 

( 0.01 ) 

13.3 ( 

: o.o2) 

1.64 ( 

; 0.034 ) 

3.14 

6.9 

4.2 

IRAS 18032-2032 

11.6 

( 0.04 ) 

1.9 ( 

: o.io ) 

6.22 ( 

: 0.663 ) 

1.87 

6.2 

1.0 

IRAS 18075-1956 

5.0 

( 0.08 ) 

1.3 ( 

: o.io ) 

3.80 ( 

; 0.463 ) 




IRAS 18162-2048 

12.4 

( 0.03 ) 

0.7 ( 

; 0.08 ) 

18.65 ( 

; 1.503 ) 

1.06 

11.7 

0.6 

IRAS 18100-1854 

22.2 

( 0.01 ) 

5.3 ( 

: o.o2 ) 

4.20 ( 

; o.no) 




IRAS 18174-1612 

9.5 

( 0.02 ) 

2.6 ( 

: o.o5 ) 

3.70 ( 

; 0.183 ) 

6.23 

1.5 

0.4 

IRAS 18317-0757 

4.4 

( 0.03 ) 

0.9 ( 

: o.o7 ) 

4.69 ( 

; 0.366 ) 

0.86 

5.1 

1.1 

IRAS 18316-0602 

32.4 

( 0.01 ) 

2.8 ( 

: o.o5 ) 

11.60 ( 

; 0.625 ) 




IRAS 18403-0417 

25.1 

( 0.01 ) 

5.8 ( 

: o.o3 ) 

4.36 ( 

; 0.125 ) 




IRAS 18479-0005 

18.1 

( 0.02 ) 

2.1 ( 

: o.io ) 

8.61 ( 

; 0.838 ) 




IRAS 18507+0110 

82.0 

( 0.01 ) 

16.0 ( 

: o.o2 ) 

5.13 ( 

; 0.106 ) 

8.86 

9.3 

1.8 

IRAS 19095+0930 

12.7 

( 0.03 ) 

1.9 ( 

: o.o7 ) 

6.70 ( 

; 0.518 ) 

1.38 

9.2 

1.4 

IRAS 19110+1045 

15.9 

( 0.02 ) 

0.9 ( 

: o.i3 ) 

17.39 ( 

; 2.230 ) 

0.49 

32.1 

1.8 

IRAS20081+3122 

14.1 

( 0.02 ) 

7.1 ( 

: o.o3 ) 

1.99 ( 

; 0.069 ) 




IRAS20255+3712 

7.9 

( 0.05 ) 

1.4 ( 

: o.o7 ) 

5.52 ( 

; 0.464 ) 

1.27 

6.3 

1.1 

IRAS20178+4046 

2.3 

( 0.00 ) 

0.8 ( 

: o.o9 ) 

2.82 ( 

; 0.258 ) 

0.66 

3.5 

1.2 

IRAS20350+4126 

4.5 

( 0.05 ) 

1.2 ( 

: o.o7 ) 

3.77 ( 

; 0.309 ) 

1.01 

4.4 

1.2 

IRAS22543+6145 

26.4 

( 0.02 ) 

4.0 ( 

: o.o8 ) 

6.56 ( 

; 0.533 ) 

3.70 

7.1 

1.1 

IRAS23138+5945 

8.3 

( 0.05 ) 

1.0 ( 

: o.i3 ) 

8.31 ( 

; i-i84 ) 




IRAS23133+6050 

15.9 

( 0.03 ) 

1.7 ( 

: 0.12 ) 

9.39 ( 

; i-i5i) 

0.80 

19.9 

2.1 


Table 15. Mean abundance ratio for each phases. 


Evolutionary Phases 


Mean HCN/HNC abundance ratio 

IRDC 

Quiescent core 

0.97(±0.10) 


Active core 

2.65(±0.88) 

HMPO 


4.17(±1.03) 

UCHII 


8.96(±3.32) 










Table 16. 


Conditions for chemical modeling 


Evolutionary Phases 

n 

[cm~ 3 ] 

T d 

[K] 

T gas 

[K] 

Chemical Timescale 

[yr] 

IRDC 


27 

27 

10000 

HMPO 

10 5 

50 

50 

60000 

UCHII 


100 

100 

50000 






